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■  ...TRACT  1 

__A  geometrical  shock  wave  transmission  model  is  used  to  determine  water  surface 
peak  spall  velocities  above  single-charge  cratering  detonations  in  a  medium  over- 
lain  by  a  thin- to-moderately- thick  layer  of  water.  These  velocities,  in  conjunc¬ 
tion  with  en  acoustic  simple  source  treatment,  predict  the  water- shock- induced 
airblast  from  such  events.  Required  input  data  include  the  acoustic  impedance  and 
shock  wave  transmission  properties  of  the  cratering  medium  in  the  inelastic  and 
elastic  regions.  The  method  is  intended  to  predict  far-field  airblast  in  the  "low" 
overpressure  regime,  for  which  case  the  acoustic  simple  source  treatment  is  valid,. 
This  regime  is  of  direct  interest  in  determining  airblast  damage  from  typical  ex 
plosive  cratering  events.  The  method  is  expected  to  be  less  accurate  for  extremely 
high  overpressure  levels  and  supersonic  peak  spall  velocities,  which  require  non¬ 
acoustic  analysis. 

t>A  computer  code  called  "MULTIMEDIA"  is  included  with  this  report  as  an  aid  to 
numerical  prediction  calculations.  Predictions  made  using  the  code  are  compared 
with  the  limited  airblast  data  available  for  underwater  cratering  tests.  Results 
indicate  that  the  basic  assumptions  of  the  method  are  correct.  Further  improve¬ 
ments  in  the  prediction  scheme  and  input  parameters  will  be  possible  as  more  data 
become  available. 
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Preface 

The  U.  S.  Army  Engineer  Waterways  Experiment  Station  (USAEWES)  Explosive 
Excavation  Research  Laboratory  (EERL)  was  the  USAEWES  Explosive  Excavation 
Research  Office  (EERO)  prior  to  21  April  1972.  Prior  to  1  August  1971  the  organiza¬ 
tion  was  known  as  the  USAE  Nuclear  Cratering  Group. 

EERL  is  conducting  a  program  of  research  in  topical  areas  critical  to  the  overall 
technology  titled  "explosive  excavation."  The  prediction  of  airblast  overpressures  is 
one  of  these  topical  areas.  This  report  presents  a  shock  wave  transmission  model  that 
is  used  to  determine  peak  water  surface  spall  velocities  above  single-charge  cratering 
detonations  in  a  medium  overlain  by  a  layer  of  water.  These  spall  velocities  are  then 
used  to  predict  far -field  peak  airblast  overpressures.  The  work  was  performed  in  con¬ 
junction  with  studies  of  airblast  safety  by  the  EERL  Environmental  Effects  Division.  It 
was  funded  by  the  Office  of  the  Chief  of  Engineers  (OCE)  appropriation  96x3121,  Gen¬ 
eral  Investigations. 

The  Directors  of  EERL  during  the  preparation  of  this  report  were  LTC  R.  L. 
LaFrenz  and  LTC  R.  R.  Mills. 

The  author  wishes  to  thank  Mr.  E.  J.  Leahy,  Chief,  Environmental  Effects  Divi¬ 
sion,  for  his  assistance  and  suggestions  throughout  the  airblast  investigation  program. 
Appreciation  is  also  expressed  to  MAJ  R.  H.  Gates  for  providing  data  and  useful  dis¬ 
cussions  concerning  media  characteristics.  Mr.  B.  B.  Redpath,  Mr.  W.  C.  Day,  and 
CPT  D.  E.  Burton  provided  additional  data  concerning  properties  of  saturated  media. 

Dr.  A.  Sakurai  of  the  U.  S.  Army  Engineer  Waterways  Experiment  Station  Weapons 
Effects  Laboratory  also  presented  a  number  of  helpful  comments  and  suggested  pos¬ 
sible  areas  for  future  research. 
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Nomenclature 


W 

Dap 

Dr 

Dw 

dob 


M 

Pambient 

PSTP 

APlccal 


AP 


rock 


^P(RW> 


AP 


water 


AP, 


s 


SGZ 

S_ 


=  Undisturbed  sonic  velocity  in  air  at  sea  level  (fl/sec). 

=  Compressional  elastic  wave  velocity  in  the  cratering  medium  (ft/sec). 
=  Undisturbed  sonic  velocity  in  water  (fl/sec). 

=  Total  apparent  charge  depth,  allowing  for  refraction  at  the  rock- 
water  interface. 

-  Depth  of  burial  of  the  charge  in  the  cratering  medium  (ft). 

=  Depth  of  the  water  layer  above  the  cratering  medium  (ft). 

=  Scaled  depth  of  burial  of  a  charge  in  s  one-layer  cratering  medium 
(no  overlying  water  layer)  =  Depth  of  burial /W1^,  (expressed  in 


ftykt1''3). 

=  Vertical  height  above  the  water  surface  of  an  elevated  far-field 
point  at  which  AP  predictions  are  desired. 

=  IV’ach  number  of  the  air  shock  at  SGZ. 

=  Ambient  air  pressure  at  SGZ,  just  above  the  water  surface. 

=  Ambient  air  pressure  at  standard  sea  level  conditions  =  1013  mbar 
=  14.88  psi. 

=  Peak  local  air-blast  overpressure  (excess  above  ambient  air  pres¬ 
sure)  at  the  rising  water  surface  in  the  near  field  (R  /(D  +  D  ) 

w  r  w 

<  1  to  2). 

--  Peak  radial  stress  in  the  rod"  layer  (excess  above  ambient  pres¬ 
sure);  AProck  is  a  function  of  scaled  range  from  the  shot  point. 

=  Peak  airblast  overpressure  (excess  above  ambient  air  pressure)  at 
the  water  surface  level  in  the  far  field  (R^/fD^  +  Dw)  >  2). 

=  Peak  shock  pressure  (excess  above  ambient  hydrostatic  pressure) 
in  the  water  layer;  calculated  for  locations  just  below  the  water 
surface. 

=  Peak  local  or  near-field  airblast  overpressure  (excess  above  ambient 
air  pressure)  at  SGZ. 

=  Exponent  describing  the  dependence  of  rock  surface  vertical  spall 
velocity  on  inverse  dimensionless  slant  range  (Drfir)  from  shot 
point.  The  quantity  "'p”  applies  to  free  reck  surface  velocities  (no 
overlying  water  layer). 

=  Scaled  range  from  shot  point  in  the  rock  layer  (ft/kt*^)- 

=  Water  surface  range:  distance  from  SGZ,  measured  along  the  water- 
air  interface  (ft). 

*  Surface  ground  zero;  the  point  at  the  water-air  interface  (or  rock- 
air  interface)  vertically  above  the  shot  point. 

*  Slant  range  from  shot  point  to  a  location  on  the  rock-water  inter¬ 
face  (ft). 
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U 


r 


^radial 

Uw 

V(rock) 


Vspall  at  SGZ 
=  VQ(rock) 

V{  water) 


VQ(  water) 
W 


a 


7 

a 

elevation 


0w 

^aCa 

PrCr 


Slant  range  from  a  point  on  the  rock-water  interface  to  a  noint  «n 
the  water-air  interface  along  a  path  of  shock  propagation  (ft). 
Press,--e  transmission  factor  at  the  rock-water  interface;  ratio  of 
the  pressure  on  the  water  side  of  the  interface  to  the  pressure  on  the 
rock  side  of  the  interface  at  a  given  location. 

Velocity  of  propagation  of  the  peak  of  the  radial  stress  wave  in  the 
inelastic  region  of  the  cratering  medium  (ft/sec;  Ur  is  less  than  the 
compressional  elastic  wave  velocity  for  most  compressible  media 
of  interest  for  excavation  applications. 

Peak  radial  free-field  particle  velocity  in  the  cratering  medium 
near  an  explosion;  rSLCttal  is  a  function  °*  scaled  range  from  the 
shot  point. 

Local  shock  front  propagation  velocity  in  water  (ft /sec);  U  =  Cw  for 
peak  pressures  less  than  20,000  psi;  Uw  >  Cw  for  very  high  pres¬ 
sures. 

Peak  vertical  spall  velocity  of  a  free  surface  with  no  overlying  water 
layer.  This  velocity  is  not  necessarily  at  SGZ  (interface  vertically 
above  the  shot  point)  but  may  be  located  at  any  distance  from  SGZ. 


Peak  vertical  spall  velocity  at  SGZ  for  a  free  surface  (that  is,  for  a 
one-layer  medium  with  no  overlying  water  layer). 

Peak  vertical  spall  velocity  of  the  water  surface  layer  (as  a  function 
of  Rw). 

Peak  vertical  spall  velocity  of  the  water  surface  at  or  close  to  SGZ. 
Yield  of  a  detonation  (kt  or  tons). 


Uniquely  determined  constant  (for  a  given  detonation)  which  relates 
peak  airblast  overpressure  to  the  dimensionless  range  in  Montan’s 
airblast  prediction  method. 

Ratio  of  specific  heats  =  1.4  for  an  ideal  gas. 

Measured  elevation  angle  to  an  elevated  location;  angle  between  the 
vertical  direction  and  a  vector  from  SGZ  to  the  elevated  location. 
Angle  between  the  direction  of  shock  propagation  and  the  rock-water 
interface  on  the  rock  side  of  ihc  interface  (degress). 

Angle  between  the  direction  of  shock  propagation  and  the  rock-water 
interface  on  the  water  side  of  the  interface  (degrees). 

Acoustic  impedance  of  air  under  sea  level  conditions  =0.0184 1  (psi- 
sec/ft). 

Acoustic  impedance  of  the  cratering  medium  in  the  elastic  propaga¬ 
tion  regime  (psi-sec/rt). 

Impedance  of  the  cratering  medium  in  the  ineiastic  regime,  for 
peak  radial  stress  wave  propagation  velocity  Ur  (psi-secrfl). 
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P  C 
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Subscripts 

a 

ap 

r 

R„ 

•» 

w 


=  Acoustic  impedance  of  water  (psi-sec/ft). 


=  Air;  parameter  or  quantity  for  the  air  layer. 

=  apparent  charge  (depth),  aMowing  for  refraction. 

=  rock  or  cratering  medium;  parameter  or  quantity  for  the  rock  layer. 
=  scaled  range  from  shot  point  in  the  rock  layer  (ft/kt*^);  equivalent 
to  true  range  i>r  for  an  assumed  yield  of  1  kt. 

=  water;  parameter  or  quantity  for  the  water  layer. 


Conversion  Factors 

British  units  of  measurement  used  in  this  report  can  be  converted  to  metric  units 
as  follows: 


Multiply 

By 

To  obtain 

inches 

2.54 

centimeters 

feet 

0.30.48 

meters 

cubic  feet 

0.02832 

cubic  meters 

cubic  yards 

0.764555 

cubic  meters 

pounds 

0.4535924 

kilograms 

pounds  per  square  inch 

0.00689476 

meganewtons  per  square  meter 

pounds  per  cubic  foot 

16.02 

kilograms  per  cubic  meter 

Fahrenheit  degrees 

5/9 

Celsius  or  Kelvin  degrees3 

foot-pounds 

0.138255 

meter-kilograms 

^o  obtain  Celsius  (C)  temperauire  readings  from  Fahrenheit  (F)  readings,  use  the 
following  formula:  C  =  (5/9)  (F  -  32).  To  obtain  Kelvin  (K)  readings,  use:  K  =  (5/9) 
(F  -  32)  ♦  273  J5. 
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ESTIMATING  WATER-SHOCK-INDUCED  AIRBLAST 
FROM  DETONATIONS  IN  A  MEDIUM  OVERLAIN 

WITH  WATER 


Section  1.  Principles 

The  method  developed  in  this  report  predicts  water-shock- induced  airblast  from 
detonations  in  cratering  media  overlain  by  shallow  layers  of  water.  The  method  utilizes 
a  simplified  model  of  shock  transmission  through  the  cratering  and  water  layers  in 
which  all  pressure  pulses  are  assumed  to  propagate  according  to  the  laws  of  regular 
transmission  and  reflection.  The  pressure  pulses  propagating  through  the  water  reach 
the  water-air  interface  and  launch  a  surface  spall  layer  according  .  -  f-e  free-surface 
law  (all  incident  pulses  are  perfectly  reflected  as  tension  waves,  uv-i-lmg  the  spall 
launch  velocity).  The  spalled  surface  induces  an  airblast  pulse  in  the  overlying  air. 
Far-field  overpressures  near  the  water  surface  level  are  predicted  by  using  the  cal¬ 
culated  vertical  spall  velocities.  Basic  assumptions  implicit  in  the  method  are  dis¬ 
cussed  below. 

Section  2.  Assumptions 

A  spherically  divergent  detonation  shock  wave  propagates  through  a  two-layered 
homogeneous  quiescent  medium  (the  bottom  layer  is  referred  to  as  "rock"  for  purposes 
of  convenience;  the  second  layer  is  water).  Transmission  and  reflection  of  shock 
waves  at  the  rock- water  and  water-air  interfaces  are  described  by  the  acoustic  laws 
for  ideal  discontinuous  waves  of  infinitesimal  amplitude.  '  The  acoustic  laws  of 
regular  transmission  and  reflection  generally  hold  for  most  cases  of  interest.  In  water, 
for  example,  shock  waves  closely  follow  acoustic  behavior  at  peak  pressures  up  to 
10,000  psi  when  the  wave's  direction  of  propagation  intersects  an  interface  at  angles 
greater  than  60°  (Ref.  1).  The  peak  radial  stress  functions  within  the  cratering  medium 
are  assumed  to  be  known.  These  peak  stresses  are  In  turn  used  to  calculate  the  peak 
pressures  transmitted  across  the  rock- water  interface  into  the  water  layer.  Elastic 
precursors  in  the  rock  layer  are  ignored.  The  peak  stress  wave  is  indeed  a  very 
strongly  dominant  effect  at  ranges  close  to  explosions  in  almost  all  cratering  media  of 
interest.  Next,  peak  pressures  t:ansmitted  through  the  water  layer  to  the  water-air 

if  precursor  pulses  should  prove  dominant  in  launching  the  water  surface  spall  for 
any  case,  the  treatment  given  below  could  still  be  used.  However,  it  would  be  necessary 
(l)  to  determine  the  peak  excess  pressure  function  of  the  precursor,  (2)  to  substitute 
that  function  for  the  usual  peak  radial  stress  function  in  the  rock  layer,  &Prock  ^  Sr>, 
and  (3)  to  substitute  the  elastic  precursor  velocity  in  rock,  Cr,  for  the  inelastic  velocity 
in  the  rock,  Ur(see  Section  4). 

-I- 


interface  are  calculated  using  semiacoustic  transmission  equations,  which  include 

allowance  for  spherical  aberration  and  astigmatism  of  the  refracted  rays  (see  Section  5 

and  Appendix  A).  It  is  assumed  that  the  water-air  interface  acts  as  a  free  surface; 

the  water  pressure  pulses  are  perfectly  reflected  ?is  tension  waves,  and  the  top  portion 

3 

of  the  water  layer  immediately  spalls  upward  as  a  free  surface.  The  water  surface 
spall  mound,  in  turn,  acts  as  a  simple  acoustic  source  of  far-field  airblast  overpres¬ 
sures  transmitted  into  the  air  layer. The  entire  "direct  ray"  approach  to  calculating 
water  pressures  and  surface  spall  velocities  is  based  on  the  assumption  that  the  water 
layer  is  sufficiently  thin  and  the  cavitational  spall  launch  after  peak  pressure  pulse 
arrival  is  sufficiently  rapid  that  direct  geometric  pulse  transmission  is  a  satisfactory 
approximation.  This  assumption  is  equivalent  to  adopting  a  local  prr~sure  approxima¬ 
tion  throughout  the  second  layer.  The  local  pressure  approximation  will  be  valid  only 
as  long  as  the  water  layer  pressure  pulses  originating  from  a  broad  area  of  the  rock- 
water  interface  are  not  combined  into  a  single  or  far-field  pressure  pulse.  Combina¬ 
tion  will  not  occur  if  consideration  is  restricted  to  water  layers  not  appreciably  thicker 
than  the  charge  burial  depth  in  the  underlying  rock  layer  and  to  the  central  part  of  the 
induced  spall  velocity  field  (surface  ranges  <  a  few  times  the  charge  burial  depth). 

Section  3.  Accuracy  and  Range  of  Applicability 

The  assumption  that  transmitted  peak  stress  acts  as  a  discontinuous  pressure 
pulse  in  launching  the  spall  is  not  strictly  correct.  However,  it  is  a  useful  approxima¬ 
tion  for  obtaining  a  straightforward  mathematical  development  capable  of  predicting 
airblast.  The  deviations  of  shock  behavior  and  waveforms  in  rock  from  this  simple 
model  are  such  that  the  method  will  overestimate  peak  pi  >.  'sures  transmitted  into  the 
water  that  are  effective  in  launching  the  spall.  Thus,  the  induced  surface  velocities 
and  airblast  will  be  overpredicted.  Likewise,  the  simple  source  assumption  discussed 
above  may  not  hold  perfectly  for  the  very  broad  velocity  fields  induced  by  charges  in 
water  and  saturated  media.  However,  breadth  of  the  shockfront  tends  to  cause  in¬ 
complete  reinforcement  of  local  airblast  pulses  from  various  parts  of  the  spalled  sur¬ 
face.  Again,  the  resulting  airblast  will  be  less  than  the  predicted  level.  This 
discussion  indicates  that,  given  proper  input  parameters,  the  method  will  safely  pre¬ 
dict  the  water-shock-induced  airblast. 

The  approach  used  in  this  report  considers  only  the  transmitted  peak  stress  pulse 
and  the  induced  peak  water  surface  spall  velocities.  Therefore,  it  predicts  the  water- 
shock-induced  peak  airblast.  Events  at  much  later  times,  including  the  mound  dis¬ 
integration  and  "blowout"  (cavity  gas  vent),  are  not  considered.  Current  evidence 
indicates  that  the  gas  »«>nt  *s  very  rapidly  suppressed  with  increasing  depth  of  burial 
for  all  chemical  explosive  events  in  water  and  in  media  of  90%  or  greater  saturation. 
The  water-shock -induced  (or  "ground-shock- induced")  airblast  component  is  strongly 
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dominant  for  events  buried  more  deeply  than  170  ft^ct1^  in  water  or  a  saturated 
medium.  It  is  expected  that  this  component  will  remain  dominant  for  detonations  in 
most  cratering  media  overlain  by  water  at  total  depths  near  or  helow  170  ft^tt1^. 
Note  that  the  theory  applied  here  is  not  appropriate  for  very  shallow  detonations, 
particularly  for  those  shallower  than  about  80  ft^ct1^.  These  events  give  rise  to  an 
intense  supersoflic  water  plume  that  penetrates  the  spall  dome  at  early  time  and  to 
other  venting  effects  that  invalidate  a  simple  spall  analysis.  Shallow  events  are  dis- 
cussed  more  extensively  elsewhere,  but  consideration  here  will  be  restricted  to 
subsonic  water  interface  velocities. 


Section  4.  Input  Data 

Application  of  a  two-layer  shock  wave  transmission  model  requires  knowledge  of 
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the  peak  oressures  in  both  layers.  The  acoustic  characteristics  of  each  medium  ’ 

are  also  needed  in  crder  to  determine  shock  wave  transmission  behavior  at  the  inter¬ 
face.  Necessary  parameters  include  the  acoustic  impedances  (pC),^  which  are  known 
for  virtually  all  media.  The  fundamental  parameters  are  listed  in  Tables  1  through  4. 
Tables  1  and  2  give  the  acoustic  impedance,  pwCw,  and  the  undisturbed  sonic  velocity, 
Cw,  in  water  (as  functions  of  temperature  and  salinity).  Table  3  summarizes  a  useful 
scheme  for  classifying  excavation  media.  Table  s  lists  the  acoustic  impedances  and 
shock  yave  transmission  properties  (see  below)  of  several  excavation  media  that  may 
be  encountered  in  underwater  applications.  For  purposes  of  this  report,  the  term 
"rock"  and  the  subscripts  "rock"  and  "r"  will  be  used  for  all  cratering  media  with 
acoustic  impedances  greater  than  that  of  water. 


Table  1 .  Acoustic  impedance  of  water  at  atmospheric  pressure  for  several 
temperatures  and  salinities.7 


Salinity, 

par*s  per  thousand  (ppt)a ' 

Temperature  (°C) 

0 

TO 

ZIT 

30 

Acoustic  impedance. 

pwCw  <Psi-8ec/ft> 

0 

61.98 

63.97 

65.43 

66.49 

10 

63.08 

64.98 

66.40 

67.42 

20 

64.19 

66.04 

67.42 

63.39 

30 

65.29 

67.11 

63.43 

69.38 

35 

65.87 

67.68 

68.96 

69.85 

aSeawater  salinity  =  35  ppt. 


"References  12  through  16. 

*All  impedances  in  this  report  are  expressed  in  units  of  psi-sec/ft,  following  standard 
acoustic  usage. 
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Table  2.  Sound  velocity  of  water  at  atm^Sph&ric  pressure  for  several 
temperatures  and  salinities.^ 


Salinity,  _ Temperature  (aC) _ 

parts  per  thousand  (ppt)a  D  ”  ZT5  3TT 

Sot^nd  velocity,  Cw  (ft /sec) 


0 

4601.1 

47^8,8 

4865.8 

4956.5 

10 

4643.7 

4701,7 

4901.0 

4987.9 

20 

4687.3 

48016 

4937.3 

5020.6 

30 

4732.1 

4808,7 

4974.7 

5054.3 

35 

47  54.9 

4809.7 

4993.8 

5071.5 

aSeawater  salinity 

=  35  ppt. 

Table  3. 


Media  classification  f'OreJcrplosive  excavation. 
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I.  Media 

A.  Common  Excavation 

1.  Soil 

2.  Clay  shale 

B.  Rock  Excavation  (requires  drilling  ^nd bHasting  to  excavate) 

1 .  Weak  rock  (less  than  4000-pSi  u*eonf~ined  compressive  strength) 

2.  Intermediate-strength  rock  <400®  to  L.  6,0OO-psi  unconfined  compressive 

strength) 

3.  High-strength  rock  (greater  tha*n  i^OOO-psi  unconfined  compressive 

strength) 

II.  Lithology  (rock  excavation)  or  Soil  Clg^-sificsation  (common  excavation)  and 

and  Geologic  Structure 

III.  Degree  of  Saturation  or  Water  Content 

A.  Dry  (less  than  50%  saturated  or  le^-s  tha_n  10%  water  content— common — or 

less  than  3%  water  content — rocI«) 

B.  Wet  (between  50  and  90%  saturate  crgrreater  than  10%  water  content- 

common — or  greater  than  3 intent — rock) 

C.  Saturated  (greater  than  90%  3aturat=ed) 


Given  the  fundamental  shock  wave  pat^xaietexrs,  one  may  easily  determine  peak 
pressures  in  water.7'^  Unfortunately,  peak  r  adial—  stresses  in  cratering  media  are 
more  difficult  to  obtain.  The  peak  stress,  at  the  rock- water  interface  is 

needed  in  order  to  calculate  the  pressure  into  the  water.  Let  Sr  be  ihe 

true  slant  range  from  shot  point  to  any  givon  location  on  the  rock-water  interface. 
Then  AP  ,  will  be  expressed  as  a  function  cfS_  -  The  following  are  five  methods 
of  estimating  AP  ^  (at  Sr). 


Table  4.  Properties  of  several  media  in  which  underground  explosive  tests  have  been  conducted. 


a.  THEORETICAL  CODE  CALCULATIONS  OF  CRATERING  DYNAMICS 

Theoretical  code  calculations  of  cratering  dynamics,®  bas^d  cn  the  known 
characteristics  of  the  cratering  medium,  can  provide  values  of  ihe  peak  radial  stress. 

I 

b.  DIRECT  MEASUREMENTS  OF  PEAK  RADIAlJ  STRESS 

Direct  measurements  of  the  peak  radial  stress  for  experiments  in  the  given 
cratering  medium  would  be  most  useful.  Unfortunately,  due  to  the  enormous  difficulty 
of  performing  such  measurements  very  close  to  an  explosion,  the  available  data  are 
incomplete  for  saturated  media. 

c.  MEASURED  FREE-FIELD  PEAK  RADIAL  PARTICLE  VELOCITIES 

Measured  free-field  peak  radial  particle  velocities  in  the  given  medium  can  be 

used  to  calculate  peak  radial  stresses.  Free-field  particle  velocities  have  been 
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measured  for  a  wide  variety  of  different  media.  *  Peak  radial  particle  velocities, 
^radial*  are  we^‘^^e<*  by  a  relation  of  the  form: 

Uradial  =  K1  (true  range/W1/3)  °P  .  «> 

where  Kj  and  nr  are  constants  for  a  given  medium,  type  of  explosive,  and  propagation 

region.  Kj  and  nf  have  two  different  values  for  a  given  medium-explosive  combination; 

nr  is  relatively  large  in  the  inelastic  or  cratering  region  close  to  the  detonation.  This 

is  the  region  in  which  permanent  displacement,  crushing,  or  cracking  of  the  material 

takes  place.  At  longer  ranges,  the  stress  decreases  to  a  level  at  which  the  material 

responds  elastically.  Energy  is  not  irreversibly  deposited  in  the  medium  (although 

some  energy  degradation  continues  to  occur  due  to  natural  dissipative  processes). 

These  longer  ranges  are  referred  to  as  the  elastic  region  of  shock  propagation;  nf  is 

usually  much  smaller  in  the  elastic  region.  The  transition  between  inelastic  and 

elastic  propagation  occurs  at  (true  range/W^)  »  140  to  700  for  typical  cratering 
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media.  *  Experimental  data  show  that  4.0  >  nr>  1.0  in  the  inelastic  region,  and 
1.3  >  np  >  1.0  in  the  elastic  region  (Table  4).  Given  the  peak  radial  particle  velocity 
at  a  specified  range,  the  peak  radial  stress,  AP^^,  at  that  range  is: 

AProck  =  Uradial  (prUr}  in  1116  inelastic  region,  (2) 

*Prock  =  U radial  (prCr}  in  the  elastic  region,  (3) 

where  PrUr  =  impedance  of  the  cratering  medium  in  the  inelastic  regime,  for  peak 
radial  stress  wave  velocity  Ur;  pr =  acoustic  impedance  of  the  cratering  medium  in 
the  elastic  propagation  regime  (sonic  velocity,  Cr).  These  quantities  are  approximately 
constant  for  a  given  medium.  The  calculated  above  is  the  increase  in  stress 
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above  the  local  ambient  pressure  (ambient  pressure  is  usually  negligible  compared  to 

AProck>*  ,« 

Graphs  of  Uradjaj  and  ^procjc  as  functions  of  (true  range/W  '  )  are  given  in 

Ref.  6  for  nuclear  explosions  in  several  media.  These  experimentally  derived  graphs 

encompass  both  the  inelastic  and  elastic  regions.  A  complete  discussion  of  the  results 

is  beyond  the  scope  of  this  report,  but  values  of  or  ^Prock  f°r  niany  media 

may  be  obtained  from  the  literature. 


d.  MEASURED  PEAK  VERTICAL  SPALL  VELOCITIES 


One  means  of  approximating  uses  measured  peak  vertical  spall  velocities 

at  surface  ground  zero  (SGZ)  for  cratering  events  with  no  overlying  water  layer. 

(Peak  vertical  spall  velocities  should  not  be  confused  with  late-time  gas -acceleration- 
induced  velocities.)  Observed  peak  vertical  spall  velocities  for  experiments  in  a  given 
medium  with  a  given  type  of  explosive  are  a  function  of  scaled  depth  of  burial 
(dob  =  scaled  depth  of  burial;  for  the  one-layer  case,  with  no  overlying  water  layer, 
dob  =  rock  depth  =  Df/W1/^).  The  measured  peak  vertical  spall  velocities  of  the 
rock  at  SGZ  fit  a  relation  of  the  form; 

'Vn  at  SGZ  =  V0  (r«kl  =  K2  Wobl  •  <4> 

£ 

where  K2  and  mr  are  constants  for  a  given  medium  and  type  of  explosive.  VgpaU  at  SGZ 
relations  for  various  media  and  explosives  are  discussed  in  Ref.  5. 

Once  the  peak  vertical  spall  velocity  at  a  specified  scaled  depth  of  burial  is  known, 
it  may  be  used  to  calculate  the  peak  radial  particle  velocity  Uradijd  in  the  given  medium 
at  a  scaled  range  equal  to  the  scaled  depth  of  burial: 

U radial  {at  true  ™*e/W1/3>  =  W  [vspaU  at  SGZ  (at  the  given  dob)] 

=  1/2  [y0  (rockKat  the  given  dob]  , 

where  (true  rangc/W1^3)  =  dob  =  Dy/W1^3.  The  thus  obtained  is  substituted  into 

Eq.  (2),  giving  AProck; 


APr0ck  scaleci  range  =  dob)  =  1/2  f(prUr)  Vq  (rock]  .  (5a) 

This  method  currently  applies  to  experiments  at  cratering  depths  (but  not  to  deep  con¬ 
tained  detonations).  It  gives  Uradial  and  ^Prock  in  the  inelastic  region  (close  ranges, 
scaled  range  <  300  ft/kt1^),  and  cannot  be  used  for  longer  ranges.  The  SGZ  peak 


’’These  values  are  constants  for  the  Inelastic  region:  i.e.,  for  experiments  at  crater¬ 
ing  depths.  The  values  of  K2  and  mr  are  expected  to  change  for  the  elastic  region 
(deep  contained  events). 


vertical  spall  velocities  for  contained  experiments  are  not  always  reliable  for  deter¬ 
mining  peak  radial  particle  velocities  in  or  near  the  elastic  region.  Photographically 
measured  velocities  for  deep  experiments  tend  to  be  lower  and  attenuate  faster  than 
true  peak  particle  velocities  should,  because  the  surface  material  does  not  break  apart 
and  spall  freely,  but  responds  in  a  partially  elastic  fashion.  If  reliable  peak  surface 
velocities  for  deep  contained  experiments  become  available,  they  may  be  used  in 
Eqs.  (5a)  and  (3)  to  obtain  AProck  at  various  ranges  (substitute  PrCr  for  Ur). 

e.  MEASURED  PEAK  VERTICAL  SPALL  VELOCITY  PROFILES 


One  additional  method  for  finding  peak  pressures  in  the  rock  layer  is  available. 
This  method  is  related  to  the  procedure  described  in  the  previous  paragraph,  because 
it  utilizes  peak  vertical  spall  velocity  data  for  free  rock  surfaces  with  no  overlying 
water  layer;  it  differs  because  it  makes  use  of  measured  velocities  at  locations  other 
than  SGZ. 

Again,  let  Vq  (rock)  be  the  SGZ  peak  vertical  spall  velocity  of  a  free  rock  surface 
above  a  detonation  at  dob  =  Dr/W^.  Then,  as  was  sho\m  in  Method  (d)  above,  AP  ^ 
at  a  scaled  range  Sr/Wly^  equal  to  the  dob  is  given  by: 

AProck  [at  Sr/w3/3  =  dcb]  *  [^rV  V0  {rock)]  *  (3a) 


Measured  VQ  (rock)  values  may  be  substituted  into  this  equation  to  obtain  AP^  at 
various  ranges.  However,  the  SGZ  velocities  may  not  be  accurately  known  over  a  wide 
range  of  dob's.  In  some  cases,  the  SGZ  velocity  at  one  scaled  depth  of  burial  and  the 
manner  in  which  the  velocity  profile  decreases  moving  away  from  SGZ  are  known.  This 
information  can  also  be  used  to  derive  ^Proc£  values,  although  with  somewhat  lower 
accuracy.  First,  let  V(rock)  be  the  peak  vertical  spall  velocity  of  a  free  rock  surface 
(no  overlying  water)  at  any  location  other  than  SGZ.  Reference  5  demonstrates  that 
V(rock)  is  usually  well-fitted  by  an  empirical  relation  of  the  form: 


V(rock)  =  Vq  (rock)  J^J  , 


where  "p"  is  a  medium-dependent  exponent  that  does  not  vary  greatly  for  detonations 
at  different  scaled  depths  of  burial  in  the  same  medium.  This  equation  describes  the 
velocity  profile  of  the  spall  mound  launched  by  the  detonation. 

Next,  let  AP  ^  be  peak  excess  pressure  of  a  shock  wave  whojte  direction  of 
propagation  intersects  a  free  rock  surface  at  angle  6  .  The  spall  velocity  of  the  free 
surface  at  that  location  will  then  be: 


V(rock)  = 


or,  substituting  from  the  empirical  V(rock)  equation. 
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iProck  ’  'P  [<«rUr>  V(r“W  sBT^] 


(5b) 


Given  the  SGZ  peak  vertical  spall  velocity  Vn  (rock)  at  one  scaled  depth  of  burial 
1 A  ^ 

(dob  =  D_/W  '  )  and  the  exponent  p,  Eq.  (5b)  is  used  to  calculate  AP  ,  at  any  scaled 
r  tys  rocK 

range  (Sr/W  '  )  greater  than  the  dob.  This  equation  is  valid  within  the  range  of  in¬ 
elastic  propagation  for  a  given  medium.  If  one  knows  the  scaled  range  at  which  elastic 
propagation  takes  over,  the  peak  pressures  obtained  from  Eq.  (5b)  may  be  very 
approximately  extrapolated  into  the  elastic  region*  as  (S„/W*'^)~1  or  (Sr/W1^)-1*^. 

Equation  (5b)  must  be  considered  a  rougher  approximation  than  the  oth<r  available 
methods,  but  it  does  provide  a  useful  independent  check.  Unfortunately,  the  surface 
velocity  exponent  p  is  not  known  at  this  time  for  most  saturated  media.  Both  SGZ 
velocities  and  the  exponent  p  have  been  accurately  measured  for  saturated  clay  shale 
(VQ  =  3.8  X  105  dob-1*5  for  dob  from  140  to  200  ft/kt1^,  and  p  *  2.5).  These  values  may 
apply  to  other  saturated  weak  rock  media  as  well.  It  is  expected  that  p  will  fall  between 
2.1  and  3.0  for  most  saturated  media. 


All  of  the  above  methods  normally  give  peak  radial  particle  velocities  and/or  peak 
rock  pressures  as  a  function  of  scaled  range.  The  results  may  also  be  expressed  in 
terms  of  true  range  [true  range  =  (W^)  (scaled  range)]  when  they  are  applied  to  a 
specific  experiment.  It  is  often  convenient  to  plot  A  as  a  function  of  true  range 

or  scaled  range  for  a  given  experiment.  This  plot  will  take  the  form  of  two  straight 
lines  (Inelastic  anu  elastic  regions)  in  a  log-log  graph.  The  peak  radial  stress, 

s 

AProck»  *n  the  cratering  medium  may  be  read  directly  from  the  graph. 

Finally,  it  should  be  noted  that  all  of  the  above  methods  apply  only  for  a  cratering 
medium  identical  to  the  medium  for  the  two-layer  experiment  in  question.  If  the  crater¬ 
ing  material  under  the  water  layer  is  saturated,*  then  the  data  used  to  derive  Uracjjaj  40(5 
APpock  *nust  be  for  a  similar  saturated  medium.  This  presents  a  particular  problem  for 
Methods  (c)  and  (d),  because  the  available  data  apply  to  experiments  with  no  overlying 
water  layer.  Care  must  be  exercised  to  select  the  data  for  a  medium  as  close  to  the 
actual  experimental  conditions  as  possible  (usually,  a  completely  saturated  medium). 

*Refer  to  Appendix  B  for  further  discussion  of  shock  pressures  and  particle  velocities 
in  the  cratering  medium. 

^"Saturated”  is  defined  as  90%  or  greater  saturation  (90%  or  more  of  the  void  volume 
in  the  material  is  filled  with  water). 
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Section  5.  Predicting  Range-Dependent 
Peak  Water  Surface  Spall  Velocities 


Peek  radial  stresses  in  the  cratering  material  will  now  be  used  to  determine  the 
transmitted  peak  pressures  in  the  water  layer  and  the  resultant  spall  launch  velocities. 
Necessary  quantities  are  listed  in  Table  1  (pwCw),  Table  2  (C^),  and  Table  4  (Cr,  Ur> 

and  prUr  for  various  cratering  media).  Other  sources  may  supplement  Table  4 
if  better  data  are  available.  Record  the  appropriate  constants  for  the  test  site  in 
question  at  the  top  of  Table  5.  List  the  charge  depth  of  burial  in  the  cratering  medium, 
Dr,  and  the  depth  of  the  overlying  water  layer,  Dw,  as  well.  Figure  ;  shows  the 
experiment  configuration  and  the  nomenclature  used. 

Hie  transmitted  peak  pressures  in  the  water  layer  and  the  corresponding  water 
surface  spall  launch  velocities,  V(water),  must  be  calculated  for  several  surface 
locations.  Table  5  provides  a  convenient  format  for  tabulating  all  calculations.  First, 
the  angle  at  which  the  direction  of  shock  propagation  in  the  rock  or  cratering  layer 
intersects  the  rock-water  interface  must  be  specified.  This  angle  is  called  dr.  When 
the  shock  wave  crosses  ihe  rock- water  interface,  it  is  refracted  and  propagates  away 
from  the  interface  at  a  new  angle  0W.  The  angle  0W  may  be  determined  from  the  law 
of  refraction — Eq.  (7),  below.  After  determining  0r  and  0^,  it  is  possible  to  calculate 
the  propagation  slant  ranges  in  the  rock  layer,  S  .  and  in  the  water  layer,  S  .  With  S 
and  the  results  of  Section  4,  it  is  then  possible  to  compute  the  peak  radial  stress, 
AProck  (at  Sr>.  on  the  rock  side  of  the  rock-water  interface.  Finally,  the  pressure 
transmission  factor  T(0r)  at  the  interface  is  calculated,  and  the  pressure  transmission 
through  the  water  layer  is  determined.  The  water  pressure,  in  turn,  controls  the 
spall  launch  velocity  of  the  water  surface.  The  equations  for  calculating  all  of  these 
quantities  are  discussed  below. 

It  is  first  necessary  to  calculate  the  shock  pressures  and  the  spall  launch  velocity 
in  the  vertical  direction,  at  SGZ.  For  the  vertical  direction,  0r  =  90°  (see  the  vertical 
ray  shown  in  Fig.  1).  Because  the  vertical  ray  suffers  no  refraction  at  the  interface, 

0W  =  0r  =  90°.  The  first  line  of  Table  5  contains  spaces  for  the  0r  =  90°  calculations. 
The  last  column  has  a  space  for  Vq  (water),  the  peak  vertical  spall  velocity  at  SGZ. 

The  subsequent  lines  of  Table  5  are  used  for  recording  calculations  at  other  0r  values 
(0r  =  80°,  70°,  60°,  etc.).  A  peak  vertical  spall  velocity,  V(water),  must  be  calculated 
for  each  of  these  0r  values.  The  final  result  of  these  calculations  will  be  a  complete 
spall  velocity  profile,  encompassing  a  broad  range  of  distances  moving  away  from  SGZ, 

The  detailed  procedures  and  equations  for  calculating  the  quantities  in  Table  5 
are  as  follows: 


First,  specify  the  angle  at  which  the  direction  of  shock  propagation  in  the  rock 
layer  intersects  the  rock-water  interface,  8r.  Use  0r  to  determine  the  propagation 


slant  range,  S  ,  in  the  rock  layer: 


sin0~' 
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Calculation  of  peak  water  surface  spall  velocities  as  a  function  of  surface  range  from  SGZ 


Pa  Co  “0*0^Wl-2^ —  V-  (water)  =  SGZ  vertical  spall  velocity 

i  of  water  surface 

(at  sea  level) 


Air 

-  R  *  D  /tone  +  D  /ton 0 

SGZ 

V  (water)  =  vertical  spall 
velocity  at 

r~7  ron9cRw 


D  f  AP  (at  or  near  SGZ) 
w  J  water 


p  c  *  68.9  gyyg 

ww  *  ft 

Water!  (for  seawater 
-  at  20®C) 


Sw“Dw/sin 


.  D  /  ton  0  . 
r  r 


W/“P water  (a*  ""8*  Rw> 


COJ  6w  =  ^Cw/ CrJ  601  ®r 


AProck  ~ P^CatrangeT) 

(at  vertical  location, 
range  Df) 

S  *  D  /  sin  © 

/  r  r  r 


'Shot  point 


Fig.  1.  Shockwave  transmission  for  cratering  medium  overlain  by  water  layer  of 

small -to- moderate  thickness  (D  ~  D  ). 

w  r 


Unless  it  is  traveling  in  the  vertical  direction  (0  =  0  ~  90°),  the  shock  wave  will  be 

r  w 

refracted  upon  crossing  the  interface;  the  angle  6  at  which  the  direction  of  propagation 
enters  the  water  layer  is  given  by  the  law  cf  refraction: 


cos  0  =  cos  0 
w 


C 

w 

r  V 


If  the  propagation  range,  Sr,  is  in  the  elastic  region  for  the  given  calculation,  substi¬ 
tute  C_  for  Ur  in  Eq.  (7}  (see  Section  4).  Calculate  and  record  0  .  Input  data  for  thii 
step  may  be  obtained  from  Tables  2  and  4,  or  other  sources. 


Determine  the  propagation  slant  range,  S  ,  through  the  water  layer: 

s  .-i- 

w  sin  0_.  * 


Using  the  results  of  Section  4,  find  tne  peak  redial  stress,  APrQc^  (at  Sr),  on  the 
rock  side  of  the  rock-water  interface.  Record  AP  ^  fet  Sr)  in  the  table. 


Subject  to  the  condition  that  Cr/Cw  or  Ur/Cw  >  1.0.  From  Huy  gen's  principle,  the 
law  cf  refraction  (Eq.  (7)]  applies  if  cos  0r  £  Cr/Cw  or  Ur/Cw,  which  is  always  true  if 
the  velocity  ratio  is  greater  than  l. 
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Next,  calculate  the  pressure  transmission  factor  at  the  rock-water  interface, 

T,  which  may  be  expressed  as  T(0r)  for  a  specified  set  of  input  parameteis.  The 
value  of  T(0  )  (ignoring  shear  in  the  cratering  medium  and  considering  only  the  stress 

j}C  *■ 

wave  )  is  given  by: 


2  PwCw  cos<90‘  '  «r> 

T(6r>  =  p~Vw  cosISO  -  er)-+prirr-cos(yo  '--^  * 


Table  6. 


Acoustic  impedance,  ppC  ,  for 
typical  ^underwater  cratering 
media.' 


Medium 

p  C 

Hr  r 

(psi-sec/ft) 

Silicon  sand 

144 

Coral  sand 

118 

Sea  sanda 

214 

Coral 

136-206 

Granite 

734 

Rock  (cement) 

374 

aVineyard  sound. 


where  the  angle  0W  is  obtained  from 
Eq.  (7),  above.  If  the  propagation  range 
Sr  is  in  the  elastic  region,  substitute 
pr Cr  for  pr Ur  in  Eq.  (9),  as  usual.  Typ¬ 
ical  values  of  p  U  and  p.C  are  given  in 
r  r  r  r 

Table  4,  while  values  of  PrCr  for  addi¬ 
tional  underwater  cratering  media  are 
listed  in  Table  6.  Table  1  summarizes 

p,  ,C  values, 
w  w 

Next,  calculate  the  water  surface 
range  from  SGZ,  at  which  the  shock 
wave  reaches  the  water  surface  {see 
Fig.  1): 


Rw  =  laiTF 


w 


tan  0 


w 


(10) 


R  obviously  approaches  zero  for  6  -  SO0  (vertical  ray). 

W  T 

Determine  the  peak  pressure  in  the  water  layer  at  R  .  where  the  shock  wave 

Q  11  15  ” 

reaches  the  water-air  interface" *  9  : 


APwater  (at  surface  ran3e  Rw)  ~ 


T<at0r)AProck(atSr) 


e  >* 

* 

" 1 . 

ian(90°  -  6  ) 
r 

tan(90  -  0  )[  1  +tan2{90  -  0  )] 
D  t  D  w  -  w 

► 

• 

^  r  w  tan(90  -  0  )[  1  +  tan2 (90  -  0  )J  J 

- 

1.13 


.  iila) 


- 5 - 

Shear  pulses  are  almost  always  negligible  for  completely  saturated  media,  which 
have  Poisson  ratios  close  to  1/2,  similar  to  a  perfect  fluid. 

^This  and  subsequent  equations  use  an  acoustic  approximation  and  the  local  sonic 
velocity,  C.  for  water.  However,  it  can  be  shown*  that  the  shock  front  velocity  is 
nearly  equal  to  Cw  and  that  the  acoustic  equations  apply  in  most  cases  of  interest  for 
water  shock  pressures  AP  ^er  <  20,000  psi. 

Refer  to  discussion  in  Appendix  A. 
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For  the  special  case  of  8r  =  90‘  (vertical  direction),  this  equation  reduces  to: 
APwater  (at  8urface  ran8e  Rw  =  °>  = 

r  °  i1,13 

Tbt  er  =  90-!  iProck  (a.  Dr>  [c,r  T  IC>r)fOjJ  • 


(lib) 


If  the  propagation  range  (Dr  in  this  case)  is  in  the  elastic  region,  substitute  Cr  for  Ur 
inEq.  (lib). 

Finally,  the  peak  vertical  spall  velocity,  V(water),  at  surface  range  R  is  given 
by: 

2  APwater  (at  surface  range  R^) 


V(water)  = 


TT~ 

Hw  w 


sin  8  . 

w 


(12a) 


For  the  special  case  of  6r  -  90°,  this  equation  reduces  to: 

2  Ap .  (at  surface  range  R,„  =  0) 

V0(water)  = - ^ -  (12b) 

where  V^(water)  is  the  peak  vertical  spall  velocity  at  SGZ. 

Perform  the  above  calculations  for  each  listed  in  Table  5;  record  the  cal- 

r 

culated  R^  and  V(water)  values  in  each  case. 

Section  6.  Predicting  Water-Shock-Induced 
Airblast  Overpressures 

The  peak  vertical  spall  velocity  profile  is  a  vertical  cross  section  of  the  detona¬ 
tion  velocity  mound  taken  through  SGZ.  Given  the  velocities  calculated  in  the  previous 
section,  it  is  possible  to  plot  the  spall  velocity  profile  and  to  predict  far-field  water- 
shock-induced  airblast  overpressures.  The  far-field  overpressures  at  surface  level 

are  assumed  to  attenuate  acoustically,  as  R  *  The  approximate  value  of  the  source 

w  4  5. 

strength  constant  in  the  acoustic  equation  is  obtained  by  Montan's  method,  *  which 
estimates  the  relative  source  strength  by  using  the  SGZ  spall  velocity  Vq  and  the 
integral  of  the  spall  velocity  profile.  It  is  first  necessary  to  integrate  the  peak  spall 
velocity  profile  graphically.  On  linear  graph  paper,  plot  the  peak  surface  spall 
velocities,  V(water),  determined  in  Section  5,  as  a  function  of  surface  range,  R^  from 
SGZ;  V(waier)  is  usually  expressed  in  ft/sec,  and  R^  in  feet,  while  SGZ  is  located  at 
Rw  =  0.0.  The  highest  velocity  is,  of  course,  VQ(water)  at  SGZ.  Flot  all  other  V(water) 
points  symmetrically  on  both  sides  of  SGZ  (in  other  words,  the  range  of  integration  for 
Rw  extends  from  -=o  to  +®;  the  gradually  decreasing  velocity  profile  extends  sym¬ 
metrically  away  from  SGZ  on  both  sides).  Draw  a  smooth  curve  through  the  plotted 
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V(water)  points.  This  curve  is  the  surface  velocity  profile.  Integrate  graphically  to 
obtain  the  integral: 

V(water)  dRw  (expressed  in 

ft  /sec) 


In  reality,  the  extreme  outer  wings  of  the  velocity  profile  have  very  little  effect  on  the 
integration  and  may  be  ignored.  Therefore,  the  integral  is  taken  over  some  reason¬ 
able  interval,  -R1  to  +R*  ,  such  that  the  important  central  part  of  the  velocity  field  is 

w  vv 

included: 


+R' 

/  " 
-Rw 


V( water)  dR 


w 


CO 

i 


V(water)  dRw  . 


Substitute  the  value  of  the  velocity  profile  integral  into  Eq.  (13)  and  calculate  "or,"  which 
is  a  dimensionless  constant  for  a  given  detonation: 


a  = 


2j r  (£>r  +  £>w)  V0  (water)  J  M 


/+» 
00 


V(water)  dR 


w 


(13) 


Next,  calculate  the  local  or  near-field  airblast  overpressure  atSGZ,  ApQ: 

Ap0  ■  W  Mv0 


(14) 


where 

Pambient  =  ambient  air  pressure  at  SGZ  for  the  given  experiment  (may  be 
determined  from  meteorological  data  or  from  the  altitude  of  the  experiment  location). 

PSTP  =  amb*en*  a*r  pressure  at  standard  sea-level  conditions;  Pqrpp  =  103  3  mbar 
=  14.68  psi. 

p  C  =  acoustic  impedance  of  air  under  standard  sea-leve)  conditions;  p_C 
=  0,01841  psi-sec/ft. 

VQ  (water)  =  peak  vertical  spall  velocity  at  SGZ  (in  ft/sec); 

M  =  Mach  number  of  the  air  shock  at  SGZ — Eq.  (15);  M  »  1.0  for  small  surface 
spall  velocities. 

If  the  above  units  are  used  for  PQCa  and  Vq  (water),  then  APq  will  be  in  psi. 

The  Mach  number  of  the  air  shock  at  SGZ  i3  determined  from  the  equation: 


M2 


Vq  (water) 


M  -  1  =  0 


(15) 


The  integral  is  considered  to  be  positive  on  both  sides  of  Rw  =  0;  in  other  words,  the 
total  integral  is  twice  the  value  of  the  profile  integral  from  Rw  =  0  to  Rw  =  «. 
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i 


*1 


«5 
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where  y  =  1  .4  (for  air  or  any  ideal  gas)  and  Ca  is  the  undisturbed  sonic  velocity  in  air 
at  sea  level;  C.  -  1087  ft/sec.  Technical^ ,  the  true  local  sonic  velocity,  C  (local), 

cl  cl 

should  be  used  instead,  but  C  (local)  ~  1087  ft/sec  for  most  cases  of  interest. 

a 

Equation  (15)  is  rather  tedious  to  solve  by  hand;  the  solution  is  plotted  as  a  func¬ 
tion  of  <V0/Ca)  in  Ref.  5. 

The  constants  APq  and  a  calculated  above  are  used  to  predict  the  far-field  air- 
blast  overpressure,  AP,  as  a  function  of  surface  range,  R^: 


AP  (at  range  Rw>  =  APQ  a 


AP  and  APq  in  psi. 


(16)* 


1 

I 


* 

f 


4 


I 


Equation  (16)  applies  only  to  overpressures  at  or  just  above  the  water  surface  level, 
and  only  to  far-field  overpressures  at  ranges  beyond  the  region  of  appreciable  water 
surface  spall  motion.  The  far-field  for  broad  velocity  profiles  of  the  type  in  question 
here5  is  normally  defined  by: 


R 

+WD  " =  2 
r  w 


(17) 


Predict  the  overpressures,  AP,  at  several  far-field  ranges,  R^.  Plot  the  resultant 
points  on  log-log  graph  paper  and  connect  them  with  a  straight  line.  The  line  has  a 
slope  of  rJ'°  and  may  be  used  to  predict  airblast  at  any  far-field  range  closer  than 
the  region  of  appreciable  meteorological  effects  (several  miles  from  SGZ  for  large- 
yield  events;  slightly  closer  for  small  yields). 

At  near-field  ranges  (inside  the  region  of  appreciable  water  surface  spall  motion), 
the  overpressure  pulses  from  various  parts  of  the  rising  water  surface  do  not  have  an 
opportunity  to  combine.  Therefore,  the  uncombined  or  local  spall-induced  overpres¬ 
sure  is  a  quantity  of  interest.  The  local  overpressure  at  any  near-field  location  is 
given  by: 

P 

APlocal  (at  range  V  -T?1— 6—  <PaCa>  MV  (water,  at  range  Rw>,  (18) 

STP 

where  M  is  a  function  of  the  local  velocity  ratio,  V(water)/C  — refer  to  Eq.  (15). 

a 

Equation  (18)  applies  only  for  locations  at  or  just  above  the  water  surface  level,  and 
only  for  near-field  ranges  (Rw/Dr  +  Dw)  <  1  to  2.  The  local  overpressure  at  SGZ,  -*Pq» 
has  already  been  calculated — Eq.  (14).  Local  overpressures  at  all  other  ranges  should 
be  correctly  predicted  by  Eq.  (18).  However,  some  combination  and  reinforcement  of 
overpressure  pulses  can  occur  even  at  near-field  ranges.  For  this  reason,  APq  may 
be  used  as  a  safe-sided  estimate  of  the  true  peak  overpressure  anywhere  in  the  near¬ 
field. 

— y 

"More  correctly,  one  should  use  apparent  charge  depth.  Pap,  (allowing  for  refraction 
of  the  wavefront)  rather  than  the  true  depth  (D  +  D^)  in  this  equation.  However,  terms 
involving  D^p  cancel  between  Eqs.  (13)  and  (lb),  allowing  the  use  of  any  desired  depth 
normalization  (see  Appendix  C). 


* 

4 
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The  above  airblast  predictions 
apply  to  locations  at  or  a  few  feet  above 
the  water  surface.  An  approximate  tech¬ 
nique  developed  by  Montan  is  sometimes 
used  to  predict  far-field  overpressures  at 
greater  elevations.  Consider  a  location 
at  height,  h,  above  the  water  surface 
(Fig.  2).  Drop  a  perpendicular  from  this 
location  to  the  water  surface.  The 
perpendicular  intersects  the  water  at  a 
surface  range  Rw  from  SGZ.  Draw  a 
line  from  SGZ  to  the  elevated  location. 
The  length  of  this  line  is  the  true  range: 


True  range 


The  line  intersects  the  perpendicular  at 
an  angle  "©elevation*"  given  by 


cos  e 


elevation 


ion‘hfel+hj112-  <20! 


Elevated  location 


“elevation 


SGZ 

coseelevatian=h/^h2  +  Rw2)1/2 
=  h/ (range  to  elevated  location) 


Fig.  2.  Elevated  airblast  (above  water 
surface  level). 


The  peal:  overpressure  at  the  elevated  location,  Apeievate(j»  is  then4: 


elevated 


[at  r”*>  (r2w + h2)I/2]  •  ^0  * 


^elevation' 


Equation  (21)  is  valid  only  for  far-field  ranges  +  h^1^/(Dr  +  Dw)J  >  2,  and  only  for 
angles  ©election5*  15°.  Even  within  this  region,  the  equation  is  not  rigorously  correct. 
Its  use  is  justified  primarily  by  the  fact  that  it  produces  accurate  overpressure  pre¬ 
dictions  for  intermediate  far-field  ranges  and  angles  ©eievation  not  extremely  close  to 
the  vertical.5  Rigorous  derivation  of  the  far-field  elevated  overpressures  would  require 
a  complete  knowledge  of  the  velocity-time  history  of  the  entire  water  surface  spall. 

Such  information  is  not  available  for  the  complex  multipie-layer  cases  under  discussion 
here. 


Section  7.  Sample  Problem 

Perhaps  the  most  interesting  case  for  applied  airblast  predictions  is  saturated 
high-strength  rock  overlain  by  water.  The  acoustic  impedance  and  shock  transmission 
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properties  of  the  cratering  medium  differ  greatly  from  those  of  water  in  this  case. 

Appropriate  input  data  for  most  saturated  high-strength  media  are  somewhat  lacking. 

However,  the  shock  transmission  properties  of  saturated  andesite  have  been  extensively 

studied.  The  available  data  derive  from  nuclear  experiments,  but  may  be  provisionally 

applied  to  all  types  of  events  until  better  information  is  obtained. 

Observations  of  nuclear  tests  in  saturated  andesite  have  shown  that  the  free-field 

6 

peak  radial  particle  velocities  fit  empirical  relations  of  the  form  - 
Inelastic  region,  Rg  <  580  ft/kt1/3:  Uradial  =  7.81  X  105  Rg1*8, 

Elastic  region,  Rg  >  580  ft/kt1^:  Uradial  =  4.70  X  103  R"1,0. 

where  R  is  the  scaled  range  from  shot  point  in  ft/kt1^3. 

The  acoustic  impedances  and  the  peak  radial  stress  wave  and  sonic  wave  velocities 
for  saturated  andesite  are  listed  in  Table  4: 

Inelastic  region:  PrUr  =  300  psi-sec/ft,  Ur  =  9500  fi/sec. 

Elastic  region:  PrCr  =  395  psi-sec/ft,  Cr  =  12,500  ft/sec. 

Hie  above  information  is  substituted  into  Eqs.  (2)  and  (3),  giving  the  peak 
stresses  in  the  rock  layer: 

Inelastic  region:  Rg  <  580  ftyktlj/3:  AProd  =  2.342  X  108  Rg1-8, 

Elastic  region:  Rg  >  580  ftjkt1''3:  *Prock  =  1.856  X  106  Rg1,0, 


where  AProc^  is  the  peak  radial  stress  (psi).  The  acoustic  impedance  and  sonic 
velocity  in  the  overlying  water  layer  are  assumed  to  be: 

Water:  PWCW  =  68.96  psi-sec/ft,  Cw  =  4994  ftyisec  (see  Tables  1  and  2,  seawater 
at  20°C). 


These  data  may  be  used  to  calculate  peak  water  surface  spall  velocities  by  the 
method  discussed  in  Section  5.  The  AProc^  equations  are  expressed  in  units  of  scaled 
range  (f^ct*^3);  therefore,  it  is  convenient  to  use  scaled  distance  units  throughout  the 
prediction  calculations.  This  procedure  does  not  limit  the  applicability  of  the  results, 
since  the  predicted  Ap^  and  a  values  will  apply  to  all  detonations  at  the  same  scaled 
rock  and  water  depth.  Only  the  distance  units  in  the  final  airblast  prediction — see 
Eqs.  (16)  or  (21) — will  change  for  events  that  differ  from  1  kt  in  yield. 

As  a  sample  proliem,  a  typical  prediction  cnlculation  Las  been  performed  for 

rock  depth  D  =  200  t*/k t^,  and  overlying  water  depth  D  =  60  ftykt^3.  The  results  of 

r  w  1/3 

this  calculation,  extending  to  a  water  surface  range  of  771  ftykt  '  from  SGZ,  are  shown 

in  Table  7.  Note  that  all  dimensions  listed  in  the  body  of  the  table  are  scaled  (ftyct1^). 

The  predicted  peal,  vertical  spall  velocities  of  the  water  surface  are  given  in  the  last 

column  of  the  tabli.  The  (maximum)  velocity  at  SGZ  :s  Vq  (water)  =  155.3  fysec.  The 

velocity  profile  may  be  integrated,  as  discussed  in  Section  6,  to  give  a.  Velocities 

from  Table  7  may  be  plotted  on  linear  graph  paper  as  a  function  of  surface  range  Rw 

from  SGZ  (profile  plotted  symmetrically,  on  both  sides  of  SGZ).  The  area  under  the 

resultant  curve  may  then  be  integrated  graphically  to  some  reasonable  range  from 
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Table  7.  Sample  problem.  Calculation  of  peak  water  surface  spall  velocities  as  a  function  of  surface  range  from  SGZ. 
prCr  =  395  psi-sec/ft  Cr  =  12,500  ft/sec  Pwcw  =  68‘96  pal-sec/Et  Cw  =  4994  ft/sec 

prUr  =  300  psi-sec/ft  Ur  =  9,500  ft/sec  Dw  =  60  ft/kt1^ 

D  =  200  ft/kt1/S 


SGZ  (to  Rw  =  771  on  each  side  of  SGZ  in  this  case).  The  profile  integral  is  used 

to  determine  a . 

rw=+771  ft/kt1/3 

V(water)  dR  =  75,230  scaled  ft2/sec  (by  graphical  integration) 

w=-771  f^tt1^ 

a  =  +T5wr V0Twater)  f  V(water)  dRw  =  ERZBOjVsSjaT  {75'230)  =  °-297* 


The  local  airblast  overpressure  APq  at  SGZ  may  also  be  calculated — Eq.  (14). 
First,  assume  an  ambient  air  pressure  of  PatT,hl-pnt  =  1013  mbar  (sea  level).  The 
acoustic  impedance  of  air  is  paC„  -  0.01841  psi-sec/ft.  The  Mach  number  M,  as 
determined  from  Eq.  (15),  is  M  =  1.089.  The  local  overpressure  APq  is  then 
APq  =  (PaTT,h}pnt/PSTP)  (Paca)MV0  (water)  =  3.11  psi.  This  result  may  be  corrected  to 
any  other  ambient  air  pressure  by  substituting  the  appropriate  value  for  P ambient* 

It  is  now  possible  to  apply  Eq.  (16)  and  to  predict  the  far-field  airblast  over¬ 
pressures  AP: 


AP  =AP0o 


=  3.11  (0.297) 
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where  AP  is  expressed  in  psi,  and  Rw  is  in  scaled  units,  ft/kt1/3.  Note  that  the  above 
equation  applies  to  locations  at  or  near  water  surface  level,  and  for  far-field  ranges 
(Rw/Dr  +  Dw)  >  2,  or  Rw  j>  520  ft^ct1^.  Overpressures  at  elevated  locations  ,r;ay  be 
predicted  using  Eq.  (21),  subject  to  the  limitations  stated  in  Section  6. 

The  scaled  predictions  above  may  be  easily  adapted  for  any  desired  explosive 
yield.  As  an  example,  perform  an  airblast  prediction  for  a  5-ton  (0.005-kt)  equivalent 
yield  chemical  explosive  charge  under  the  conditions  of  tbe  sample  problem  (change 
scaled  dimensions  to  true  dimensions): 


Dp  =  200  flykt1^3  =  34  ft, 
Dw  =  60  ftyct1/3  =  10  ft. 


and,  from  Eq.  (16), 


AP  =  AP0« 


=  3.11  (0.297) 


AP  =  2-1-  where  R  is  in  true  (instead  of  scaled)  feet, 

w 

The  prediction  is  valid  near  water  surface  level  and  for  ranges  R^  £  88  ft. 
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Note  that  all  airblast  predictions  in  this  report  are  based  on  "acoustic"  or  R~*‘® 
overpressure  attenuation.  Effects  such  as  deviation  from  acoustic  behavior  (strong 
shocks),  nonsymmetric  shock  propagation,  and  atmospheric  temperature  gradient  may 
cause  the  attenuation  rate  to  vary.  The  wide  changes  in  observed  attenuation  rates  for 
similar  experiments  indicate  that  local  effects  (rather  than  nonacoustic  shock  effects) 
are  responsible  for  deviations  from  R-1'®  attenuation.  Temperature  gradient  is  the 
most  important  local  factor  influencing  far-field  airblast.  At  locations  near  the  ocean, 
the  lower  layers  of  air  typically  show  a  strong  temperature  gradient;  such  conditions 
cause  sound  to  be  refracted  upwards,  away  from  the  surface  level.  Thus,  the  over¬ 
pressure  pulse  is  expected  to  attenuate  somewhat  more  rapidly  than  R*^*®  at  surface 
level.  If  a  strong  gradient  prevails,  the  overpressure  predictions  may  be  slightly 

improved  by  attaching  a  line  of  slope  R-1"2  to  the  normal  R*1*^  line  at  "long"  ranges, 

w  r  W  -12 

beyond  the  region  of  close-in  acoustic  attenuation  [very  roughtly,  attach  an  *  line 
at  ranges  greater  than  Rw  =  10  to  20  times  (Dr  +  Dw)J.  This  procedure  gives  an 
approximation  of  the  faster  attenuation  and  lower  overpressures  at  ranges  beyond  the 
close-in  region. 

It  should  also  be  noted  that  the  opposite  atmospheric  effect  of  a  temperature 

j$c 

inversion  can  cause  refraction  of  sound  waves  toward  ground  level,  with  resultant 
focusing.  Under  these  conditions,  overpressures  may  attenuate  more  slowly  than 
R^  *®,  producing  increased  airblast  at  long  ranges  (=  a  few  miles  from  SGZ)  where 
the  focusing  is  most  likely  to  occur.  Other  atmospheric  effects  may  focus  sound  waves 
at  very  long  ranges,  but  these  effects  become  important  only  for  large  multi-kiloton 
detonations  capable  of  creating  damage  at  such  extreme  ranges 


Section  8.  Summary 

A  method  is  presented  for  estimating  water-shock -induced  airblast  from 
single-charge  cratering  detonations  in  saturated  media  overlain  with  water.  The 
method  requires  as  input  data  a  knowledge  of  the  acoustic  impedances  and  sonic 
velocities  in  the  cratering  medium  and  water  layers,  as  well  as  an  empirical  fit  to 
the  shock  propagation  parameters  in  the  cratering  medium.  Predictions  made  using 
this  technique  have  not  yet  been  adequately  compared  with  actual  airblast  data,  and 
future  observations  may  reveal  the  need  for  some  improvements.  To  verify  and  to 
extend  the  method,  measurements  of  airblast  and  time-dependent  water  surface  spall 
velocity  for  an  underwater  cratering  event  in  intermediate-  or  high-strength  rock  are 
needed. 

The  procedures  described  in  this  report  are  easily  adapted  for  computer  use. 

A  relatively  simple  code  may  be  written  to  calculate  peak  water  surface  spall  velocities 
as  a  function  of  range  R^,  over  small  intervals  in  R^  (or  over  small  increments  of  ©r). 

Temperature  increases  with  altitude  in  the  near-surface  layers. 
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The  resultant  velocity  profile  may  be  integrated  numerically  to  some  reasonable  outer 
limit  in  or  V(water),*  giving  the  value  of  a  .  Computer  calculations  provide  a 
particularly  convenient  means  of  calculating  APq  and  a  for  a  number  of  different  rock 
and  water  depth  combinations  in  a  variety  of  media.  A  computer  program  that  cal¬ 
culates  APq  and  a  for  depth  combinations  20  <  Dff  <  320  ft/kt^  and  60  <  Dr  <  320  ft^ct1^* 
is  presented  in  Appendix  B.  Some  sample  problems  solved  using  this  program  are 
discussed. 


$ 

A  reasonable  criterion  for  limiting  the  calculation  might  be  R  =  5(D  +  D  ),  or 
V(water) «  0.02  V0(water).  w  r  w 


) 

I 
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Appendix  A 

Pulse  Transmission  in  the  Water  Layer 
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Appendix  A 

Pulse  Transmission  in  the  Water  Layer 


An  approximation  that  describes  the  elastic  propagation  of  refracted  acoustic 
compression  waves  in  two-layered  media  has  been  developed  by  Hipp,  Gerjouy,  Horton, 
and  others.* This  approximation,  derived  from  fundamental  considerations  of 
geometry  and  energy  conservation,  allows  for  spherical  aberration  and  astigmatism  of 
the  wavefront  in  the  second  medium.  It  does  not  allow  for  curvature  of  the  wavefront 
or  nonacoustic  pressure  attenuation.  With  the  insertion  of  a  small  correction  for 
attenuation,  the  Hipp-Gerjouy-Horton  treatment  is  used  in  this  report.  Complete 
derivations  may  be  found  in  the  literature,  but  the  equations  are  summarized  here  for 
the  reader's  convenience. 

Assume  that  the  second  layer  consists  of  water,  and  let  AP^  be  the  peak  water 
pressure  caused  by  a  given  ray  at  the  rock- water  interface.  Then  AP^.  may  be  cal¬ 
culated  from  Eq.  (9): 


=  T(er>*Prock<atSr>' 


where 

T(0  _)  =  acoustic  transmission  factor  at  the  rock-water  interface 
r 

AProck  (at  Sr)  =  rock  pressure  at  the  rock-water  interface. 

The  refracted  compression  wave,  AP^.,  will  continue  to  propagate  through  the 
water  layer,  reaching  the  surface  at  surface  range,  R^,  from  SGZ.  The  water  pres¬ 
sure  at  the  surface,  APTO^er  (at  surface  range  R^),  is  then  given  by: 

f  D 

APwater  (at  surface  ran*e  V  =  APwi  p? 


X 


tan(90°  -  0  ) 

_ r _ 

tan(90  -  0  )[1  +  tan2(90  -  0  )) 

D  +  D  w  w 

•  r  wtan(90  -  ©r)[l  +tan3(90  -  0r>]  . 


The  A-term  is  a  correction  factor  that  allows  for  nonacoustic  pressure  attenuation  in 
the  second  layer;  A  •»  1.13  for  water.  The  above  equation  is  in  turn  used  to  derive 
Eq.  (11)  of  the  text. 


-24- 


/ 


Appendix  B 

MULTIMEDIA-Computer  Calculation 
of  Two-Layer  Airblaet  Parameters 


Appendix  B 

MULTIMEDIA— Computer  Calculation 
of  Two-Layer  Airblast  Parameters 

MULTIMEDIA  is  a  computer  program  (code)  that  utilizes  the  method  discussed 
in  the  foregoing  report  to  calculate  the  peak  vertical  spall  velocity  profile  of  the  water 
surface  above  detonations  in  two-layered  media  (cratering  medium  overlain  by  water), 
to  integrate  the  velocity  profile,  and  to  determine  a  and  AP^.  The  quantities  a  and 
APq  may  then  be  used  to  predict  peak  airblast  overpressures,  AP,  as  described  in 
Section  6. 

The  code  calculates  necessary  shock  wave  propagation  parameters  for  waves 
intersecting  the  rock- water  interface  at  1°  intervals,  starting  with  9r  =  89°  (©r  =  89°, 

88°,  87°,  86°,  . . .  ).  The  calculated  shock  wave  parameters  at  6r  =  89°  are  used  to 
approximate  the  water  pressure  and  the  peak  vertical  water  surface  spall  velocity 
VQ  (water)  at  SGZ.  This  approximation  is  adopted  as  a  matter  of  simplicity  and 
convenience  in  the  computer  code;  it  has  no  significant  effect  on  the  numerical  accuracy 
of  the  results.  The  velocity,  Vq,  is  in  turn  used  to  obtain  the  peak  local  airblast  over¬ 
pressure,  APq  (in  psi),  at  SGZ.  The  calculated  shock  wave  parameters  for  other  0f 
angles  give  the  water  surface  spall  velocities,  V(  water),  as  a  function  of  surface  range, 
R^,  from  SGZ.  The  calculations  are  terminated  when  the  peak  vertical  water  surface 
spall  velocity  decreases  below  -jr^th  of  Vq  (or  when  0r  =  4°,  whichever  occurs  first). 

The  velocity  profile  of  the  water  surface  is  integrated  to  the  termination  point  sym¬ 
metrically  on  both  sides  of  SGZ  (the  velocity  V(water)  at  0r  =  39°  is  again  used  as  an 
estimate  of  VQ  (water)  in  the  integration].  The  integral  of  the  velocity  profile  is  then 
used  to  calculate  the  value  of  a — see  Eq.  (13)  of  the  text. 

$ 

A  complete  list  of  the  code  is  included  at  the  end  of  this  appendix. 

INSTRUCTIONS 

In  order  to  use  the  code,  it  is  necessary  to  know  the  acoustic  impedance  (in 
psi-sec/ft)  and  sonic  velocity  (in  ft/sec)  in  the  water  layer,  acoustic  impedance  and 
sonic  velocity  for  both  the  inelastic  and  elastic  regions  of  the  cratering  medium,  the 
scaled  range  of  transition  between  the  inelastic  and  elastic  regions,  the  peak  radial 
particle  velocities  caused  by  the  shock  wave  in  the  inelastic  and  elastic  regions,  and 
the  ambient  air  pressure  at  the  experiment  location  (in  millibars).  Determination  of 
these  input  quantities  is  discussed  in  the  text  of  this  repor*.  Read  the  data  into  the 
computer  in  the  following  format: 

First  data  card:  All  eighty  columns:  in  alphabetical  and  numerical  characters, 
enter  the  name  of  the  medium,  ambient  air  pressure,  and  any  other  desired  information. 

*The  MULTIMEDIA  code  is  written  in  FORTRAN  II  and  runs  approximately  3  min 
on  a  CDC  6600  computer. 
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Second  card:  Enter  the  acoustic  impedance  (psi-sec/ft)  of  water  in  the  first  five 
columns  (format  F5.0),  and  the  sonic  velocity  (ft/sec)  In  the  water  in  the  second  five 
columns  (F5.0). 

Third  card:  Enter  the  acoustic  impedance  (psi-sec/ft)  in  the  inelastic  region  of 
the  cratering  medium  in  the  first  five  columns  (F5.0),  the  stress  wave  velocity  (ft/sec) 
in  the  inelastic  region  of  the  cratering  medium  in  the  next  five  columns  (F5.0),  the 
scaled  range  of  transition  between  the  inelastic  and  elastic  regions  (in  ftyfct1^3)  in  the 
next  five  columns  (F5.0),  the  acoustic  impedance  in  the  elastic  region  in  the  next  five 
columns  (F5.0),  and  the  sonic  velocity  in  the  elastic  region  in  the  next  five  columns 
(F5.0). 

Fourth  card:  The  peak  radial  particle  velocities  in  the  cratering  medium  due  to 
the  shock  wave  are  fitted  by  two  connected  (log-log  straight  line)  equations  of  the  form: 

B1 

Inelastic  region:  Vparticje  (in  ft/sec)  =  (scaled  range)  , 

Bo 

Elastic  region:  Vparticje  (in  ft/sec)  =  Ag  (scaled  range)  , 

where  "scaled  range"  is  the  distance  from  the  charge  center,  expressed  in  ft/kt1^. 

In  the  first  nine  columns  of  the  data  card  enter  A^  (format  E9.0),  in  the  next  five 
columns  enter  B,  (format  F5.0),  in  the  next  five  columns  enter  the  scaled  range  cf 
transition  between  the  inelastic  and  elastic  regions  (ft^ct1^)  (format  F5.0),  in  the  next 
nine  columns  enter  Ag  (format  E9.0),  and  in  the  next  five  columns  enter  Bg  (format 
F5.0). 

Fifth  card:  Enter  the  ambient  air  pressure  at  the  experiment  location  (millibars) 
in  the  first  eight  columns  (format  F8.2).  The  ambient  pressure  at  sea  level  is  usually 
assumed  to  be  1013.0  mbar,  The  ambient  pressure  at  5000  ft  MSL.  is  approximately 
844.0  mbar. 


CODE  OUTPUT 


The  code  output  occurs  in  the  following  format:  At  the  top  of  the  first  output  page, 
all  of  the  identifying  information  punched  on  the  first  input  card  is  printed  out.  This 
information  includes  the  medium  name,  ambient  air  pressure,  and  anything  else 
appearing  on  the  first  data  card.  Beneath  the  identifying  information,  calculated  shock 
wave  parameters  for  the  first  combination  of  water  and  rock  depth  are  printed  out.  The 
first  depth  combination  consists  of  a  water  depth  Dw  =  20  ft/kt1^  and  a  rock  depth 
Dr  =  60  ft/kt1^.  The  parameters  are  given  for  1°  intervals  in  angle  of  intersection  with 
the.  rock  interface,  beginning  at  Gr  =  89°  (83°,  87°,  . . . ).  In  the  first  column,  the  angle 
at  which  the  direction  of  shock  wave  propagation  intersects  the  rock  side  of  the  rock- 
water  interface,  0r,  is  given;  in  the  next  column,  the  angle  at  which  the  (refracted) 
shock  wave  enters  the  water  layer,  0^  is  listed;  in  the  third  column,  the  surface 

1  to  " 

range  R  (ftykt  from  SGZ)  at  which  the  shock  wave  reaches  the  water  surface  is 
w 
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listed;  in  the  fourth  column,  the  peak  rock  pressure  at  the  rock- water  interface  is 
given.  The  next  two  columns  list  the  transmission  factor  into  the  water  layer  T(0r) 
and  the  peak  water  pressure  at  the  water  surface,  (at  Rw)  (in  psi).  In  the 

seventh  column,  the  peak  vertical  spall  velocity  of  the  water  surface  (in  ft/sec)  is 
listed;  the  final  column  gives  the  cumulative  velocity  profile  integral  (in  scaled  ft 2 /sec) 
to  this  surface  range  on  both  sides  of  SGZ.  These  shock  wave  parameters  are  cal¬ 
culated  and  printed  until  the  peak  vertical  spall  velocity  decreases  below  l/50th  of  the 

V„  at  SGZ.  The  calculation  is  then  concluded,  and  the  final  results  are  printed  two 
0  1/3 

lines  below:  first,  the  scaled  water  depth  and  rock  depth  (ft^t  '  )  are  printed  out, 

followed  by  the  Mach  number,  M,  of  the  air  shock  at  SGZ,  then  the  approximate  peak 

vertical  spall  velocity,  Vq,  of  the  water  surface  at  SGZ  (ft/sec),  the  peak  local  airblast 

overpressure  APq  at  SGZ  (psi),  and  finally  the  airblast  alpha  value  a  (dimensionless). 

The  local  overpressure  at  SGZ,  AP^,  and  the  dimensionless  constant,  a,  may  be  used 

to  predict  water-shock-induced  airblast  for  the  given  combination  of  water  depth  and 

rock  depth,  as  discussed  in  Section  6  of  the  paper. 

Next,  the  calculated  shock  wave  parameters  for  the  second  combination  of  water 
depth  and  rock  depth  are  printed  out.  The  rock  depth  is  increased  by  20  ft/kt1'^  over 
the  previous  value  in  this  case.  The  output  format  is  exactly  the  same  as  above: 
first,  the  calculated  shock  wave  parameters  are  given  for  1-deg  intervals  in  angle  of 
intersection  with  the  rock  interface,  starting  at  0r  =  89°.  The  calculation  is  again 
concluded  when  the  vertical  spall  velocity  decreases  below  l/50th  of  the  VQ  at  SGZ. 

Final  results  are  printed  two  lines  below  (the  given  scaled  water  depth  and  rock  depth, 
Mach  number  M,  peak  vertical  spall  velocity  VQ  at  SGZ,  peak  local  overpressure 

Aprt  at  SGZ,  and  the  dimensionless  constant  a). 

"  1/3 

The  rock  depth  is  next  increased  by  another  20  ftykt  '  ,  and  results  for  this  depth 

are  printed  out.  In  this  fashion,  the  program  continues  to  increase  the  rock  depth  at 
20  ft/kt1^  intervals  until  a  depth  of  320  ft^t1^3  is  reached.  Complete  output  results  are 
printed  for  each  depth.  After  a  depth  of  320  ft/kt1^  is  reached,  the  program  returns 
the  rock  depth  to  its  original  (smallest)  value  of  60  ft^ct^  and  increases  the  water 
depth  by  20  ft^ct1^3.  Complete  calculations  are  performed  for  this  new  depth  combina¬ 
tion  and  output  results  are  printed  in  the  usual  format.  The  rock  depth  is  again  varied 
at  20  ft/kt1^3  intervals  up  to  320  ft/kt1^3,  and  complete  results  are  printed  for  each  depth. 
Afterwards,  the  water  depth  is  increased  by  yet  another  20  ft/kt1^3,  and  complete  data 
are  calculated  and  printed  out  for  each  rock  depth.  This  process  continues  until  the 
final  water  depth  of  320  ft^ct^  is  reached.  The  last  set  of  calculations  and  the  last 
set  of  printed  data  output  are  for  the  deepest  combination  of  rock  and  water  depth  (rock 
depth  =  320  ft/kt1^3  and  water  depth  =  320  ft/kt1^3 ). 

For  certain  input  parameters  and  certain  combinations  or  rock  and  water  depth, 
the  physical  discontinuity  between  inelastic  and  elastic  propagation  regions  in  the  rock 
may  cause  the  change  in  surface  range,  dR^  to  become  negative  over  one  angle 
interval.  This  angle  interval  always  occurs  for  the  ray  at  the  slant  range,  Sr,  where 
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the  elastic  region  first  begins.  When  a  negative  change  of  R  is  encountered,  the 
program  prints  a  warning  immediately  below  the  data  for  the  angle  d.c  where  this  change 
takes  place.  Such  a  warning  will  occur  no  more  than  once  in  the  tabulated  data  for  a 
given  depth  combination.  The  program  uses  an  absolute  value  of  dRw  to  calculate  the 
velocity  profile  integral  contribution  over  this  interval.  Actually,  the  physical  dis¬ 
continuity  has  no  significant  numerical  effect  on  the  final  results;  thus,  it  may  be 
ignored  for  the  purposes  of  interest  here. 

Finally,  after  all  the  depth  combinations  have  been  fully  calculated,  the 
important  prediction  parameters  APq  and  a  are  summarized  in  a  convenient  table  at 
the  end  of  the  program  output.  Horizontally  at  the  top  of  this  table,  all  the  water  depths 
(Dw  =  20  to  320  ftykt1^)  are  listed.  Vertically  at  the  left  side,  the  rock  depths 
(Dr  =  60  to  320  ftykt1^)  are  given.  Within  the  table,  APq  and  a  values  for  each  depth 
combination  are  compiled.  The  APq  (in  psi)  for  each  depth  combination  is  printed 
directly  above  the  a  value.  To  read  a  value  from  this  table,  enter  at  Dw  on  the  top 
list,  proceed  down  to  Dr  (listed  at  the  left),  and  read  the  desired  APq  (top  number)  and 
a  (bottom  number). 

To  perform  overpressure  predictions  for  any  given  case,  select  values  of  APq 
and  a  for  the  set  of  data  closest  to  the  desired  combination  of  rock  depth  and  water 
depth.  Use  these  values  to  predict  water-shock-induced  airblast.  In  some  cases,  it 
may  be  necessary  to  interpolate  linearly  for  APQ  and  a  between  tabulated  values  at 
rock  depths  or  water  depths  differing  by  20  ft^ct1'  . 

SPECIAL  CASES;  MODIFICATIONS  TO  THE  CODE  AND  PREDICTION  PROCEDURES 

The  output  of  the  MULTIMEDIA  code  or  the  cede  itself  may  be  modified  to  allow 
for  certain  variations  in  the  specified  conditions.  Foi  example,  the  predicted  SGZ 
local  airblast  overpressures  APq  are  based  on  the  assumed  input  value  of  ambient  air 
pressure,  Pambient  ^nPut^-  (See  the  discussion  of  input  data  cards,  above.)  After 
running  the  program,  it  may  become  necessary  to  convert  the  calculated  APq  values 
to  new  APq  values  for  a  different  ambient  pressure,  Pambient  (new).  Any  APQ  value 
calculated  by  the  program  may  be  approximately  converted  to  a  r.ew  APQ  value  accord¬ 
ing  to  the  following  formula: 

P  ,  .  ,  (new) 

NewAPC  =  P*™' "  (lnpu£7*P0  <calculated>* 
ambient  e 

The  a  values  are  not  changed  by  a  new  value  of  Pambienf 

Also,  the  water  and  rock  depths  and  depth  intervals  for  which  quantities  are 
calculated  may  be  changed  from  the  values  specified  in  the  program.  Statements  28 
and  32  control  the  water  and  rock  depths,  respectively.  These  are  ”D#)"  statements 
that  tell  the  program  to  perform  calculations  from  a  specified  minimum  value  to  a 
specified  maximum  value  (for  example:  "D (f)  27  JJ  =  3,  16";  in  this  case,  3  is  the 
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minimum  value  and  1 6  is  the  maximum  value).  The  first  or  smallest  depth  for  which 
calculations  are  performed  is  twenty  times  the  minimum  value;  the  last  or  greatest 
depth  for  which  calculations  are  performed  is  twenty  times  the  maximum  value.  To 
change  the  minimum  depth.,  modify  the  first  (minimum)  value  in  the  appropriate  state¬ 
ment;  to  change  the  maximum  depth,  modify  the  second  (maximum)  value  in  the 
appropriate  statement.  Calculations  will  still  be  performed  at  20  ftykt1^3  intervals 

jje 

between  the  minimum  and  maximum  depths,  as  usual.  If  the  minimum  rock  depth  is 
changed  during  this  step  (Statement  32),  Statement  122  must  also  be  changed  to  agree 
with  the  new  minimum  depth.  Statement  122  says  "D$?04  JK  =  3,  JKMAX."  The 
number  "3"  must  be  changed  to  l/20th  of  the  new  minimum  depth. t 

Note  that  calculations  performed  by  the  code  for  very'  shallow  depths  may  be 
inaccurate  due  to  the  high  water  surface  velocities  (Vq  approaches  or  exceeds  the 
local  sonic  velocity  in  air).  In  addition,  shallow  explosions  may  produce  very  strong 
supersonic  plume  and  gas-vent  effects,  and  these  effects  may  become  dominant  in 

t 0-1 fi  * 

causing  far-field  airblast.  Thus,  the  code  output  cannotrbe  used  to  perform  air- 

blast  predictions  for  very  shallow  scaled  depths.  It  should  ai^o  be  noted  that  code 
calculations  at  extremely  deep  scaled  depths  cf  burst  (total  depth  »  several  hundred 
ftykt1^)  are  meaningless  for  acoustic  airblast  predictions,  because  perfect  acoustic 
reinforcement  between  pulses  from  various  parts  of  the  rising  surface  no  longer  occurs. 
The  true  far-field  overpressures  fall  below  predicted  acoustic  values,  and  local  over¬ 
pressures  may  become  important  for  some  safety  considerations.  This  situation 
occurs,  however,  only  for  very  deep  events. 

The  sonic  velocity  in  air  assumed  for  MULTIMEDIA  is  =  1087.0  ft/sec.  This 
velocity  is  usually  correct  within  about  57c  for  near  sea  level  conditions.  The  sonic 
velocity  does  vary  with  temperature,  and  a  more  accurate  vat  ue  may  he  inserted  if 
desired.  Statement  80  controls  the  sonic  velocity:  "RATIC5  =  V$/1087.0."  To  change 
the  assumed  velocity,  substitute  the  correct  value  in  the  denominator  of  this  equation. 
The  velocity  should  be  expressed  in  ft /sec. 

In  certain  cases,  there  may  be  no  clearly  defined  transition  range  between  the 
inelastic  and  elastic  shock  propagation  regimes,  or  the  medium  properties  may  not  be 
sufficiently  well  known  to  permit  the  establishment  of  two  separate  regimes.  Under 
these  conditions,  the  assumed  shock  velocities,  acoustic  impedances,  and/or  shock 
transmission  properties  may  be  identical  for  both  regions.  However,  it  is  still  neces¬ 
sary  to  repeat  the  required  input  data  twice  on  the  appropriate  input  cards.  Choose 
any  arbitrary  value  for  the  transition  range,  and  insert  this  value  on  input  cards  three 
and  four.  List  the  required  input  parameters  twice  on  these  two  cards,  for  the 
arbitrary  "inelastic"  and  "elastic"  regions  (even  though  these  data  are  the  same  for 
both  regions).  Note  that  this  procedure  is  intended  primarily  for  special  cases,  and 

sjc 

Statements  29  and  33  control  the  water  and  rock  depth  intervals  between  successive 
calculations  (factor  of  20).  These  intervals  may  also  be  modified  if  desired. 

^Numbers  specified  in  "D<?"  statements  must  always  be  integers  (no  decimal). 
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may  not  produce  accurate  results  if  the  input  data  are  used  for  regions  outside  their 
range  of  applicability  (i.e„  low-amplitude  long-range  shock  transmission  parameters 
and  long-range  "acoustic"  medium  properties  applied  to  the  close-range  inelastic 
region).  Normally,  the  transition  range  used  by  the  program  should  correspond  rather 
closely  to  a  physically  real  transition  in  the  shock  propagation  properties  of  the  medium. 

Sometimes  the  known  characteristics  of  shock  propagation  in  che  cratering 
medium  may  include  the  peak  shock  stress,  AProck,  instead  of  the  peak  particle 
velocity,  Vparticle.  The  stress  is  usually  fitted  by  two  empirical  equations  of  the 
form  ''Ap  ^  =  K  (scaled  range)®,"  one  each  for  the  inelastic  and  elastic  propagation 
regions.  The  code  requires  peak  particle  velocities  rather  than  stresses.  However, 
the  empirical  stress  equations  may  be  easily  converted  to  particle  velocity  equations: 
divide  the  constant  K  by  the  appropriate  acoustic  impedance  PrUr  or  PrCr  in  each 
region  (inelastic  and  elastic).  The  new  equations  thus  obtained  will  be  the  correct 
particle  velocity  equations  to  use  in  each  region.  Two  separate  equations  must,  of 
course,  be  calculated  (for  the  inelastic  and  elastic  regions).  The  attenuation  rates 
"B"  in  each  region,  and  the  scaled  range  of  transition  between  inelastic  and  elastic 
regions,  are  identical  for  both  stress  and  particle  velocity. 

One  additional  modification  may  be  applied  to  the  MULTIMEDIA  code,  either 
before  or  after  it  is  run.  The  integral  solutions  for  a  sometimes  converge  rather 
slowly,  giving  rise  to  large  values  of  a.  These  values  may  become  a  little  unrealistic, 
particularly  for  events  at  great  scaled  depths,  because  of  the  very  broad  linear  extent 
over  which  the  surface  velocity  field  is  integrated.  The  code  itself  may  be  modified  to 
terminate  the  integration  sooner.  Statement  102  specifies  the  termination  criterion. 

The  calculation  may  be  terminated  closer  to  SGZ  (at  a  higher  peak  velocity)  by  de¬ 
creasing  the  constant  in  this  statement.  For  example:  the  original  statement  reads 
"VTEST  =  V(p-  50.0  X  VEL."  If  it  is  desired  to  terminate  the  integral  at  a  velocity 
value  of  l/30th  VQ  (instead  of  l/50th  Vq),  the  statement  should  be  modified  to  read 
"VTEST  =  V<J)-  30.0  X  VEL."  Alternatively,  a  completely  different  termination  criterion 
may  be  established,  using  "VTEST"  as  the  decision  parameter  (see  Statement  103). 
Positive  "VTEST"  terminates  the  integral. 

Termination  of  the  velocity  integral  may  also  be  accomplished  after  the  program 
has  been  run.  The  data  output  for  each  calculation  includes  a  complete  list  of  the 
calculated  parameters  at  1°  intervals,  from  0  =  89°  to  the  termination  velocity.  The 
last  two  columns  of  this  list  give  the  Vertical  Velocity  (of  the  water  surface)  and  the 
Cumulative  Velocity  Profile  Integral  (up  to  the  listed  velocity  value).  The  final  value 
of  a,  printed  below  the  list,  is  determined  for  the  entire  integral  t  velocity  profile. 

To  terminate  the  integral  at  a  higher  velocity,  select  any  desired  termination 


& 

The  particle  velocities  given  by  the  two  equations  at  the  transition  range  should  be 
the  same;  if  they  are  not,  adiust  the  constants  slightly  as  necessary. 
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velocity  from  the  list.  Read  the  Cumulative  Velocity  Profile  Integral  at  this  velocity 
(last  column).  Calculate  a  new  value  of  a  for  this  termination  point: 


New  a  = 


1 

2  jtVq  (scaled  rock  depth  +  scaled  water  depth) 


where 


X  (CUM.  VEL.  PROFILE  INTEGRAL), 


Vq  =  vertical  velocity  of  the  water  surface  at  SGZ  (ft/sec); 
Scaled  depths  in  ft^ct1^. 


SAMPLE  PROBLEMS 


MULTIMEDIA  was  first  tested  by  application  to  certain  simplified  problems  with 
known  answers,  and  by  comparative  manual  calculations  for  selected  problems  that 
were  more  complex.  The  numerical  results  are  sufficiently  accurate  for  current 
predictive  purposes.  The  code  was  next  applied  to  several  sample  media,  using 
estimated  shock  wave  transmission  parameters  for  these  media.  The  preliminary 
trial  cases  included  saturated  andesite  (a  high-strength  rock  with  known  shock  wave 
transmission  characteristics  for  nuclear  detonations),  saturated  clay  shale  (a  weak 
rock  with  known  shock  wave  transmission  characteristics  for  nitromethane  detonations), 
saturated  coral  (a  very  weak  medium  that  fails  in  crushing,  with  estimated  shock  wave 
transmission  characteristics  for  aluminized  ammonium  nitrate  slurry — AANS — detona¬ 
tions),  and  saturated  basalt  (a  high-strength  rock  with  estimated  shock  wave  transmis¬ 
sion  characteristics  for  TNT).  The  shock  wave  transmission  characteristics  are  most 
accurately  applicable  to  the  given  explosive-medium  combination  in  each  case,  but  the 
results  will  serve  temporarily  for  estimating  airblast  from  all  sorts  of  explosives  in 
the  listed  media.  Better  approximations  to  the  shock  wave  and  medium  characteristics 
may  be  used  in  the  future,  when  more  data  from  explosive  tests  and  output  from 
SOC-TENSOR  analyses  become  available.  Input  data  for  the  various  sample  media 
are  shown  on  a  standard  computer  coding  form  (Fig.  Bl).  The  first  set  of  data  is  for 
saturated  andesite  (nuclear  events),  the  second  set  is  for  saturated  clay  shale 
(nitromethane),  the  third  set  is  for  saturated  coral  (AANS),  and  the  fourth  set  is  for 
saturated  basalt  (TNT).  Input  data  for  the  saturated  basalt  case  must  be  considered 
rough  estimates  at  this  time;  thus,  predictions  made  using  these  data  may  be  somewhat 
less  accurate  than  the  other  cases. 

As  an  example  of  the  input  format,  examine  the  list  for  saturated  clay  shale 
(second  set  of  data).  The  first  card  gives  the  medium  name  and  identification.  The 


One-  and  two-dimensional  hydrodynamic  codes  in  use  at  the  Lawrence  Livermore 
Laboratory.” 
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PROGRAMMER _ _  DME 


Fig.  B-l.  Standard  computer  coding  form. 


second  card  lists  the  properties  of  seawater  (at  20°C,  35-ppt  salinity  for  the  sample 
problem): 

Acoustic  impedance  =  68.96  psi-sec/ft 

Sonic  velocity  =  4994  ft/sec. 

The  third  card  lists  acoustic  impedances  and  shock  wave  velocities  in  the  cratering 
medium,  as  well  as  the  transition  range  between  inelastic  and  elastic  propagation 
regimes: 

Inelastic:  Acoustic  impedance  =  191.0  psi-sec/ft 
Velocity  =  6400  ft/sec 
Transition  range  =  140.0  ft^ct1^ 

Elastic:  Acoustic  impedance  =  191.0  psi-sec/ft 
Velocity  =  6400  ft/sec 

(These  quantities  do  not  differ  significantly  between  the  propagation  regions  for  highly 
saturated  clay  shale  at  ranges  greater  than  80  to  90  ft^ct1^.) 

The  fourth  card  gives  the  constants  in  the  empirical  particle  velocity  equations 
(nitromethane  in  saturated  clay  shale,  inelastic  and  elastic  regions),  as  well  as  the 
transition  range: 

Inelastic:  V  =  1.454  X  106  (scaled  range)-1*84 

Transition  range  =  140.0  ftykt1^ 

s  -1  42 

Elastic:  Vparticle  =  1.795  X  10°  (scaled  range)  *H 

The  final  card  specifies  the  ambient  air  pressure  in  mbar: 

Pambient  =  1013*°°  mbar 

The  above  shock  wave  and  medium  parameters  are  based  on  SOC-TENSOR  code 
calculations  for  chemical  explosive  cratering  detonations  near  optimum  depth  of  burial 

17 

in  saturated  clay  shale.  Code  predictions  provided  accurate  estimates  of  the  true 
peak  pressures  at  the  shock  front.  These  peak  pressures  were  converted  to  peak 
particle  velocities  via  division  b>  the  acoustic  impedance  pC  =  191.0.  Two 
straight  line  segments  were  fitted  through  the  peak  particle  velocities  for  use  by 
MULTIMEDIA.  These  straight  line  segments  accurately  approximate  al'.  of  the 
SOC-TENSOR  data  in  the  region  80  ft^ct^<  scaled  range  <  several  hundred  ft^ct1^. 
The  accepted  value  of  acoustic  impedance,  pC  =  191  psi-sec/ft,  is  also  known  to  be 
valid  for  scaled  ranges  >  90  ft^ct1^.  MULTIMEDIA  simply  converts  the  fitted  particle 
velocities  back  to  peak  pressures  for  use  during  the  calculations.  This  is  currently 
the  most  accurate  method  for  estimating  shock  pressures  in  the  clay  shale.  It  should 
be  noted  that  the  two  fitted  particle  velocity  line  segments  do  not  correspond  exactly 
to  inelastic  and  elastic  propagation  regimes  in  this  case;  they  are  the  best  empirical 
approximations  to  shock  properties  in  "high"  and  'low"  pressure  regimes  for  a  weak 
highly  saturated  medium.  Correlation  between  the  empirical  propagation  regimes  and 
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the  inelastic -elastic  transition  is  generally  much  closer  for  high-strength  media  or 
for  media  that  fail  in  crushing. 

Data  abstracted  from  the  code  output  for  die  four  sample  problems  are  given  in 
Tables  B1  through  B4  (sample  of  code  output  with  labels  added).  The  listed  informa¬ 
tion  includes  APq  (psi)  and  a  for  each  combination  water  and  rock  depth.  The  water 
depths  (ft^t1^)  are  listed  horizontally  at  the  top  of  the  table,  and  rock  depths  (ft/kt^) 
are  listed  vertically  on  the  left  side.  For  each  depth  combination,  the  program 
t  tabulates  a  value  of  AP^  (psi;  top  number)  and  a  value  of  a  (dimensionless;  bottom 

number).  It  is  important  to  remember  that  the  calculations  may  be  somewhat  in  error 
(due  to  the  high  velocities  involved)  and  that  gas -vent  airblast  may  become  dominant 
for  very  shallow  total  depths  (upper  left  corner  of  the  tables).  Final  results  for  each 
medium  may  be  read  from  the  appropriate  table.  As  an  example,  examine  Table  B1 
(saturated  high-strength  andesite):  at  water  depth  =  60  ft^ct^,  rock  depth  =  200  ft^ct1^, 
APq  =  3.11  psi,  a  =  0.309.  (This  value  of  a  differs  slightly  from  the  hand-calculated 
value  in  the  Section  7  sample  problem  because  the  integral  was  carried  only  to 
Rw  =  771  ft/kt1/3  in  the  hand-calculated  problem.)  The  values  of  APQ  and  a  read  from 
Tables  B1  to  B4  may  be  used  to  predict  far-field  water-shock-induced  airblast  over¬ 
pressures. 

,  Results  in  the  third  sample  table  are  applicable  to  saturated  coral,  and  may  be 

CIO 

compared  with  the  Tugboat  single-charge  airblast  measurements.  (The  tabulated 

solutions  are  based  on  estimated  shock  transmission  parameters  for  slurry  explosive 

* 

1  in  coral.  )  Tables  B5  through  B7  list  the  ranges  an’d  far-field  overpressures,  AP,  for 

five  Tugboat  single-charge  events.  Theoretical  a's  and  AP^'s  interpolated  from 
Table  B3  are  also  given  for  each  experiment.  The  observed  values  and  theoretical 
predictions  are  compared  in  Figs.  B2  through  B4.  Theoretical  lines  are  shown  as 
Rw  1‘°,  and  as  1  * 2  (extrapolated  from  Rw/(Dr  +  D^)  =  20.0).  Thefaster  R^,  over¬ 
pressure  attenuation  is  typical  of  temperature  gradient  conditions  in  sea  air,  which 
often  cause  sound  to  be  refracted  away  from  ground  level,  decreasing  the  observed 
overpressures.  Agreement  between  measurements  and  the  theoretical  predictions  is 
;  excellent  at  "close"  far-field  ranges  (Rw/(Dr  +  Dw)  =  5  to  30).  At  longer  ranges,  the 

far-field  overpressures  tend  to  attenuate  as  R^*^.  The  observed  overpressures  for 
Event  Id  fall  somewhat  below  the  predictions.  This  fact  could  indicate  a  slightly 
defective  shock  transmission  model  in  the  coral,  or  an  unusually  strong  atmospheric 
temperature  gradient  causing  rapid  attenuation  of  overpressure  at  the  time  of  the  shot. 

Saturated  coral  as  a  shock  transmitting  medium  does  not  vary  too  dramatically 
from  pure  seawater;  therefore,  the  above  results  do  not  provide  definitive  support  for 
the  method.  Airblast  measurements  for  a  single-charge  event  in  a  saturated  high- 
strength  medium  overlain  by  water  are  urgently  needed  in  order  to  make  further  im¬ 
provements  in  predictive  capability.  The  necessary  data  should  include  observations 

I  The  estimated  shock  wave  parameters  were  derived  from  theoretical  considera¬ 

tions.  Explosive  properties  for  aluminized  ammonium  nitrate  slurry  were  used  in  the 
derivation. 
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Table  Bl.  AP0  and  a  for  nuclear  explosives  in  saturated  andesite  (high-strength  rock) 
overlain  by  seawater  (Pambient  =  1013  mbar). 

Entries  in  pairs:  APq — in  psi,  listed  above;  a — dimensionless,  listed  below, 

SUMMARY  0T  DILI  A  VICO  AND  ALPHA  VALUES 


HATEP  DEPTHS’FT^KTl^?  LISTED  WRIJOffTALLY  ACRUSS  TOP 
WO:  DEPTHS'  FT/KTl-'SHISTED  VERTIC^.LV 


20.0 

i*0-0 

bO.O 

SO.O 

100.0 

120.0 

140.0 

lbO.O 

1X0.0 

200.0 

220.0 

2*10.0 

2bO.  0 

290.0 

300.0 

320.0 

bO.O 

49.52 

0.268 

37.6b 
0.2  72 

30.84 

0.255 

25.7*1 

0.243 

22.04 

0.234 

16.17 

0.22b 

lb.  62 
0.22-J 

15.11 

0.21b 

13.  b2 
0.211 

12.38 

0.208 

11.34 

0.204 

10.451 

0.201 

6.b8 

0.200 

8.01 

0.167 

8.42 

0.165 

r.ou 

0.164 

80.0 

2  $.77 
0.307 

16.66 

0.25*1 

17.17 

0.2b7 

15.01 

0.255 

13.26 

0.245 

11.61 

0.237 

10.77 

0.232 

6.82 

0.22b 

6.02 

0.221 

8.33 

0.217 

7.73 

0.214 

7.21 

0.211 

b.  75 
0.208 

b.  3*1 
0.20b 

5.67 

0.204 

5.1.5 

O.iot 

100.0 

14.26 

0.311 

12.52 

0.262 

11.12 

0.27fc 

*4.6*1  I 
0.285  | 

6.  OS 
0.255 

8.2b 

0.247 

7,bO 

0.240 

7.02 

0.23b 

b.  S3 
0.231 

b.  06 
0.22b 

5.71 

0. 222 

5.37 

0.216 

5.07 

0.21b 

4.  *6 
0.213 

*1.55 
0. 710 

**.32 

0.209 

120.0 

6.b7 

0.315 

3.8*4 

0.263 

7.63 

0.264 

7.21 

0.272 

b.b3 

0.253 

b.  14 
0.255 

5.71 

0.248 

5.33 

0.243 

5.00 

0.233 

4.70 

0.233 

4. *44 
0.226 

‘1.20 

0.227 

3.66 

0.223 

3. 76 
0.220 

5.b2 

0.2L7 

3.45 

0.215 

1*40.0 

7.0*4 

0.320 

b.*P4 

0.303 

5.43 

0.23*1 

0.44 

0.280 

5.11 

0.270 

4.77 

0.2b2 

4.48 

0.257 

4.22 

0.250 

3.63 

0.243 

3.77 

0.240 

3.55 
0. 237 

3.40 

0.224 

3.24 

0.230 

3.10 
0.  777 

2.67 

0.224 

2.9*1 

0.223 

lbO.O 

5.30 

0.323 

5.00 

0.30*4 

4.k5 

0.2*48 

4.35 

0.28b 

4.03 

0.273 

3.S4 

0.270 

S.b- 

0.2b3 

3.43 

0.257 

3. 2o 
0.235 

S.JO 

0.248 

2.6b 

0.244 

7.83 
0. 240 

2.70 

0.237 

:.56 

0.234 

2.46 

0.233 

2.36 

0.230 

l°0.0 

*4.20 

0.326 

4.01 

0.3W 

3.7b 

0.302 

3.54 

0.2*43 

3.34 

0.284 

3.17 

0.27b 

3.01 

0.270 

2.8b 

0.285 

2.73 

0.250 

2.bl 

0.255 

2.50 

0.251 

2.36 

0.247 

2.30 

0,.2*lb 

2.21 

0.245 

2.13 

0.240 

2.05 

0.237 

200.0 

3.50 

0.332 

3.30 

0.318 

3.11 

0.303 

2.  *45 
0.26*4 

2.80 

0.260 

2.bb 

0.282 

?  Si 

0.273 

2.43 

0.272 

2.32 

0.2bb 

2.23 

0.2b2 

2.14 

0.258 

2.0b 

0.25b 

1.68 

0.252 

1.61 
0. 246 

1.84 
0. 24b 

1.79 

0.24U 

270.0 

2.  *42 
0.338 

2.7? 

0.325 

2.b3 

0.314 

2.50 

0.3CW 

2.38 

0-263 

2.28 

0.261 

2.18 
0.  ISM 

2.0-4 

0.273 

2.01 

0.273 

1.63 

0.2L8 

l.8b 

0.2tb 

1.76 

0.2b3 

1.75 

0.256 

J.b7 

0.25b 

l.b2 

0.253 

1.5b 

0.230 

2*40.0 

2.48 

0.342 

2.&> 

0.330 

2.25 

0.31*1 

2.15 

0.313 

2.0b 

C.3C4 

1.67 

0.267 

1.36 

0.2*0 

1.82 

0.284 

1.73 

0.276 

i.b6 

0.277 

l.b3 

0.273 

1.58 
0. 2b6 

1.53 

0.*>b 

1.48 
0. 2b3 

1.43 

o.  :bO 

1.36 

0.257 

280.0 

2.14 

0.34b 

2.0*4 

3.J53 

1.65 

0.327 

1.87 
0. 318 

1.80 

0.310 

1.73 

0.303 

l.bb 

0.2*fc 

l.bO 

0.26S 

1.55 

0.288 

1.50 

0.284 

1.45 
7. 280 

1.40 

0.27b 

1.3b 

0.275 

1.32 
0. 2b6 

1.2S 

0.2bb 

1.24 

O.Tbb 

230.0 

1.8b 

0.354 

J.V* 

0.343 

1.7’ 

0.333 

1.85 

0.324 

1.56 

0.31b 

1.53 

0.306 

1.49 

0.3CH 

1.43 

0.300 

1.39 

0.265 

1.33 

0.261 

1.26 

0.297 

1.25 

0.283 

1.22 

0.276 

1.19 

0.27b 

1.15 

0.77b 

1.12 

0.273 

300.0 

1.84 

0.35*4 

1.56 

0.343 

1.52 

0.336 

1.4b 

0.330 

1.41 

0.322 

1.3b 

0.318 

1.32 
0. 312 

1.28 

0.307 

1.24 

0,302 

1.20 

0.268 

1.1b 

0.264 

1.13 

0.260 

1.10 

0.28b 

1.07 

0.28b 

1  .CM 
0. 735 

1.01 

0.230 

320.0 

1.4b 

0.3b4 

1.40 

0.354 

1.35 

0.344 

1.31 

0.340 

l.  2b 
0.332 

1.22 
0*  325 

1.18 

0.316 

1.15 

1  0.313 

1.12 

0.306 

1.08 

0.305 

1.05 

0.301 

1.02 

0.300 

1.00 
0. 267 

0. 67 
0.263 

0.6S 

0.260 

0.62 

0.288 
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Table  B2.  APq  and  a  for  chemical  explosive  events?  in  saturated  clay  shale  overlain  by 
seawater  <Pambient  =  1013  mbar). 


Entries  in  pairs:  APQ — in  psi,  listed  above;  a — dimensionless,  listed  below. 


suci**'  of  WL7 a  f»e»  ch va»  ufi 

kftTTt-  H5VTH 5>f?  *T1  '3,  1 1$T£D  K)*IZ3HTi*.LV  CCWSS  T9P 

wio:  nnj»  :r  >n,3,i uufl  vimcoxv 


r 


> 


;o.o 

40.0 

to.c- 

80.0 

100.0 

170.0 

140.0 

ito.o 

w.u 

r^.bc 

0.321 

V5. 

o.*c« 

12.15 

0.294 

?*.?5 

0.^0 

?:•.  is 

l^.ue 
r  itc 

ib.te 
0. 291 

*C  0 
35. V* 
6.324 

100.0 

20  95 
0.  341 

??.  I* 

o.i:- 

Z2.S3 

0.314 

lo.84 

U.307 

lb.  US 
0.304 

14.00 

0.299 

12.53 

0.29S 

11.04 

O.rni 

i\*4 

14.57 

0.3*" 

li-.ts 

0.321 

11.04 

0.317 

9.8l 

0.U2 

3.£2 

0.309 

3.00 

0.305 

1^0,0 

13.94 

0.  w* 

il.*? 
0.  *51 

10.30 
0. 34  *. 

4.10 

G.3ib 

3.14 

0.333 

?.» 
0.  «S 

8.89 

0.323 

b.14 

0.31S 

140.0 

9.72 

o.  ;*i 

H.+u 

0.-.'** 

i.">Z 

O.ilK 

t.ai 

0.382 

•5.10 

0.2?b 

5A5 

0.350 

!?.l9 

6.34b 

4.91 

0.342 

iW.O 

’.w, 

o.iei 

n.  5  \» 

b.3? 
0.7c  « 

5.7* 
0. 382 

5.27 

V.359 

«.&> 

0.334 

4.51 
0.  J‘X 

4.K> 

0.33*0 

ttrt.O 

t.yi 

0.2*4 

ft  it  1 

'.40 

O.a.'i 

4.-4. 

0.185 

4.:« 
0. 3Vi 

4.15 

0.3/r 

*.07 

0.353 

3.71 

0.344 

-\*.0 

'.Vi 

o.irc 

•5.00 
».  :V0 

4.8* 
0. 577 

4.33 

O.ifc? 

4.03 

0.<82 

3.  n, 
a.Jbo 

3.53 

0.53b 

3.32 

U.33? 

.ro.c 

4.81 

0.3*8 

4.44 

tt.il 

3.375 

1. 3J 
C. 

3.0 « 
0.3*4 

3.3* 

0.380 

3.l7 

0.S5S 

2.SS 

0.353 

i40.0 

M.  JV 
0.5*4, 

3.  ♦. 
•i. :  <3 

'7  * 

3.92 

0.371 

3.2? 
0.  £U> 

3.o3 

0.34.' 

2.97 

O.fU. 

2.72 

0.357 

.isO.C 

5..M 

0.38' 

5.  Vi 
3.381 

3..'? 
c.  57 

3 

2.9; 

C.MO 

2.7V: 
0.  V* 

2.01 

O.J*jO 

2.44 

O.SbO 

7*0.0 

5.  *5 

o.3^v 

*.  V. 

•  *  -  t 

?.*' 

o. 

/  -4 
fi.  '  4 

2.85 

0.170 

0.3V  0 

2.40 

0.  :t2 

2.2S 

i'.3J>2 

:0O.0 

?.oT, 

0. 

t».  v» 

: .  />• 
r».  <8>l 

C.  i75 

2.4? 
v.  571 

2.  52 

O.'b 

*•.  **1 
V 

2.11 

:to.o 

o’,  ktf 

*.»  < 
O.Jo1 

2.4." 

i'.Jol 

2.35 

0.5./ 

2.24 

C.272 

2.14 

0.289 

-'.CC 

0.3V5 

l.Sb 

0.3b? 

130.0 

200.0 

2*0.0 

240.0 

2b0.0 

290.0 

3W.0 

520.0 

14.7; 

13.25 

11.S9 

10.S2 

10.02  1 

1 

C.5S 

e.ot 

0.279 

0.27b 

0.274 

0.274 

0  272  1 

0. 271 

0.270 

0.2b3 

9.  So 

9.0b 

8.30 

7.hb 

7.  lu 

8.81 

8.1A 

5.81 

0.29A 

0.233 

0.28b 

0.284 

0.284 

0.282 

0.291 

0.279 

7.31 

b.73 

8.22 

5.79 

5.41 

5.07 

4.?7 

4.50 

0.304 

0.301 

0.293 

0.29b 

0.2=4 

0.294 

0.293 

0.291 

S.bb 

5.2b 

4.90 

4.53 

4.31 

4.0b 

3.34 

3.84 

0.319 

0.31b 

0.313 

0.311 

0.301 

0.30b 

0.307 

0.303 

4.47 

4.ie 

3.9? 

3.88 

a.  48 

3.29 

3.12 

2.9? 

0.341 

0.339 

0.335 

0.332 

O.230 

0.320 

0.320 

0.32b 

3.^3 

3.b9 

3.47 

3.28 

3.  U 

2.95 

2.81 

2.88 

0.342 

0.342 

0.334 

0.33b 

0.334 

0.332 

0.530 

0.128 

3.49 

3.29 

3.11 

2.45 

2.31 

2.87 

2.55 

2.44 

C.34i* 

0.343 

U.342 

0.340 

0.337 

0.3J5 

0.333 

0.331 

3.1- 

2.  ^7 

2.0? 

2.W 

2.55 

2.U4 

2.  S3 

2.24 

0.349 

0.34b 

0.343 

0.343 

0.341 

Q.3$e 

0.33b 

0.334 

2.54 

2.b9 

2. 5b 

2.45 

2.34 

:.:u 

2.15 

2.Cb 

0.352 

0.349 

0. 3t*b 

0.343 

0.243 

0.341 

0.539 

U.  337 

2.5* 

2.4b 

2.35 

2.25 

2.15 

2.C7 

1.93 

1.91 

0.3W 

0.351 

0.349 

0.34b 

0.343 

0.341 

0.342 

0.340 

2. 37 

2.2b 

2.17 

2.03 

1.99 

1.92 

1.P4 

1.7S 

0.33b 

0.333 

0.  *5i 

0.343 

0.34b 

0.344 

0.342 

0.342 

2.19 

?.09 

?.01 

;.43 

1.85 

1.70 

1.7? 

1.8-b 

0.  ite 

0.35b 

0.333 

0.350 

0.34* 

0.34b 

0.  4*4 

0.342 

2.02 

1.94 

1.97 

1.79 

1.73 

1.8? 

l.bi 

1.58 

0.  *>.0 

0.357 

0. 353 

0.352 

O.SSO 

0.4*3 

U.-?ib 

0.344 

1.93 

1.81 

1.74 

1.88 

1.82 

1.5b 

1-51 

l«4b 

0.359 

0.559 

0.357 

0.354 

0.352 

0.350 

0.348 

0.24b 

^Explosive  parameters  for  nitromethane  were  used  in  these  calculations. 
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Table  B3.  ApQ  and  a  for  chemical  explosive  events4  in  saturated  very  -weak  coral  over 
lain  by  seawater  <Pambient  =  1013  mbar). 

Entries  in  pairs:  APg — in  psi,  listed  above;  a — dimensionless,  listed  below. 


IflfcV  «f  WlTA  i*«6»  VCLUn 


HJiTIV  MHrHS>FT'hTl/3*  U  >TLD  WW>*0UTf4.L»  TUI* 


eoo  mrro:* 

rMTl 

•Ms&tb  vimcdiv 

ro.o 

'Ul.O 

bO.O 

%o.c 

100.0 

120.0 

140.0 

uo.o 

180.0 

2OJ.0 

220.0 

2*0.0 

0 

A 51.0 

$$0.0 

320.0 

10. 0 

43.bO 

0.33b 

51. >« 
O.V-f* 

24.2ii 

0.31K 

14.43 

0.30*4 

lb.  21 
0.305 

13.82 

0.300 

12.  Cl 

o.r*b 

10. 5^1 
o.  2*t3  j 

4.4t> 

t'.7C 

$.53 

0.20 

7.7b 

7.11 

0.28b 

b.*4- 

0.2rr$ 

b.*£ 
0.  2^3 

5.1b 

0.283 

5.28 

0.282 

00.0 

23.07 

0.340 

IS.  1" 
0.3« 

14.**tf 

0.321 

12.55 

0.31b 

10.81 

C.olO 

3.47 

O.JCfc 

3.41 

0.$02 

?  .53 
n.  vo 

t*.f« 

0.24? 

b.25 

0.285 

5.75 

O.^tl 

32 

0.242 

4. *45 
0. 2*) 

M.b2 

4.  43 
0.7J7 

4.07 

0.28». 

100.0 

14.53 

0.344 

12. Ct 
O.JSJ 

10.27 

0.32b 

8.02 
0.  $14 

7.6? 
0.  31b 

7.05 

0.311 

b.34 
•1. 507 

3. 7? 
•>.:04 

s.  .“i 
0.705 

4.37 

o.vo 

4.S2 

4.248 

4.21 

0.24b 

3.^1 

0.2^#* 

?.  TO 
0. 70  2 

?.4*t 

0.212 

3.74 

0.281 

120.0 

10.14 

0.344 

0.71 

0.32b 

7.bl 

0.330 

b.75 

0.324 

b.Ob 

0.313 

5.'*8 

0.31b 

5.01 

0.312 

4.b0 

o.  $04 

4.2^- 

O.Vjb 

3.45 

0.305 

i.M 

0.502 

3.4b 

0.300 

3.25 
0.  7*3 

?.u7 

0.2*1? 

2.*b*J 

0.7‘« 

2.75 

0.2<-5 

140.0 

7.5b 

0.347 

b.b3 

0.330 

5.-*2 

0.331 

3.3$ 

0.327 

4.84 

0.322 

4.43 

0.318 

4.oe 

0.31b 

3.73 

0.313 

3.5? 

0.310 

3.JS 
0.  50? 

?.<>» 
0.  !0.’ 

2.41 

•J. 

Ci!$02 

0.301 

2.47 

0.2-44 

2.35 

0.247 

IK>.0 

5.40 

0.344 

5.28 

0.541 

4.77 

0.3J*. 

4.34 

0.331 

3.*44 

0.72b 

3.bC 

0.524 

3.41 

0.120 

3.18 

0.?lb 

.*.48 
M.  $13 

?.e*> 

0.511 

2. 84 
0.310 

.*.50 
0.  103 

2.37 
0.  tXK 

2.2*5 

0.3*X 

2.14 

o.?o$ 

2.04 

0.301 

190.0 

4.7V 

0.350 

4.31 

0.343 

3.<M 

o.ir* 

3.b2 

0.333 

3.35 
0. 324 

3.11 

0.327 

2.*U 

4.  $2$ 

2.73 
*>.  *20 

»  ~ 

o:;t: 

2.42 

0.514 

2.24 

0.514 

2.17 
0.  $12 

2.0? 

0.310 

J.H? 
*».  f'fi 

1.88 

o.?*> 

1.80 

0.504 

200.0 

3.<W 

0.331 

J.bl 

0.547 

3.32 

0.341 

5.08 
0. 33b 

2.87 

0.33? 

2.*-3 

0.330 

2.52 

0.32b 

2.57 
k  32? 

2.24 
#*.  520 

7. 12 
0.518 

2.2? 

0.31? 

1.42 
0.  415 

t.Oj 
0. 31$ 

1.75 
0.3 11 

l.b7 

o.m 

l.bO 

0.208 

220.0 
3.3  2 
0.333 

*.07 
0.  34m 

2.fC. 

2.bb 

0.  w 

2.44 

0.334 

?.{H 

0.354 

2.  20 
o.$24 

2.LT3 
•».  32b 

1.47 

O..C3 

1.88 
0. 521 

1,74 

0.521 

1.70 

0.318 

l.b$ 

0.31b 

l.*>. 

0.31'* 

1.50 

0.31$ 

1.44 

0.311 

240.0 

2»C3 

0.357 

0.351 

2.4.* 

0.34b 

2.32 

0.341 

2.13 

0.333 

2.Cb 

0.3?b 

1.4b 

0.332 

1.05 

0.324 

1  7b 

o. 

l.bS 

0.524 

l.bO 
0.  $7* 

1.53 
0.  322 

1.47 
*).  320 

3.41 

0.318 

1.55 

o.$u> 

1.30 

0.314 

2b0.0 

7.4$ 

0.154 

2. T- 
o'rs'- 

2.17 

0.3*0 

2.C5 

0.3«3 

1.4$ 

0.342 

1.83 

0.338 

1.74 
0.  ?35 

l.V* 

«>.$$? 

1.58 

1.51 
0. 527 

1.44 

0.527 

1.33 
0.  $25 

1.4$ 
u.  $.*$ 

1.28 
u«  3.1 

1.23 

0.51*4 

1.14 

0.317 

280.0 

2.18 

O.JbO 

-r.CO 

0.355 

1.h5 

0.3o 

1.82 
0.  $«b 

1.7! 

0.345 

l.b4 

0.341 

t.% 

0.*$6 

1.44 
0. 535 

1.45 

O.S22 

I.  $7 
0.530 

1.31 

C.  5Vji 

1.2b 
0.  ZjA 

t.?l 

0.:5> 

1.  »? 
0.124 

1.12 

0. 322 

1.04 

0.320 

3CO.O 
1.4$ 
0.  Jb2 

uev 

0.357 

1.7. 

■X.V.7 

i.b4 

0.351 

1.55 

0.14? 

1.46 
0. 914 

l.Ul 
0. 341 

1.35 
u.  'J? 

l.?0 

0.3<3 

1.24 

0.335 

l.’JJ 
0. 53? 

1.15 
0.  3*1 

t.ll 

o.?^ 

1.07 

0.327 

1.03 

0.323 

1.00 

0.32$ 

C0.0 

1.73 

0.$b4 

0.  55 « 

1.5'. 

o. 

l.w$ 

0.35? 

1.41 

0.350 

1.35 

0.34b 

1.24 

0.343 

1.75 
o.  3m 

1.18 

o.$?« 

1.14 

0.538 

1.10 
0. 33b 

l.O> 
0. 334 

1.02 
0.  3$2 

fi 

o. 1 30 

0.45 

0.328 

0.82 
0.  ?2t> 

aExplosive  parameters  for  a  6%  aluminized  ammonium  nitrate  slurry  were  used  in 
these  calculations. 
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Table  B4.  APq  and  a  for  chemical  explosive  events3  in  saturated  basalt  overlain  by 
seawater  (Pambient  =  1013  mbar);  estimated  medium  and  shock  wave 
parameters  were  used  in  these  calculations. 

Entries  in  pairs:  APq — in  psi,  listed  above;  a — dimensionless,  listed  below. 

SUMMARV  ar'DELTrt  KO)  AND  ftLPhft  VALUES 


MATES  DEPTHSjTT  KT1/3.-  LISTED  HSSIHtffCTALLV  ACROSS  TOP 
MTO;  DEPTHS'  rr-rTl-'3/LISTED  VERTICALLY 


20.0 

40.0 

feO.Cl 

eo.o 

100.0 

120.  U 

140.0 

lbO.O 

1SO.O 

200.0 

220.0 

2**0. 0 

2bO.(5 

280.0 

300.0 

320.0 

bO.O 

Ll.bO 

o.ibl 

48.31 

0.2^4 

41.41 

0.215 

15.13 

0.202 

30.43 

0.181 

2b.  71 
0.183 

23.75 

0.177 

21.23 

0.172 

13.27 

0.1b8 

17.58 
0.  Ib4 

lb.  14 
O.lbl 

1*4. 80 
0. 153 

13.85 

0.155 

12.83 

0.153 

J2.Cb 

0.151 

11.32 

0.143 

iO.O 
25.  bM 
0.270 

22.17 

0.24b 

18.47 

0.223 

17.31 

0.215 

15.55 

0.205 

14.10 

0.13b 

17.87 

0.183 

11.83 

0.183 

10.34 

0.173 

10.1b 

0.175 

3.48 
0. 171 

8.83 
0.  Ib3 

3.35 
0.  Ib5 

7.87 
0.1  b2 

7.45 

0.1M5 

7.0b 

0.153 

100.0 

13.85 

0.27b 

12.  V 

0.254 

11.  S4 
0.253 

10.55 

0.225 

3.51 

0.215 

3.73 

0.20b 

8.1b 

0.133 

?.b2 

0.134 

7.14 

0.188 

b.?l 
0. 134 

b.33 
0.  173 

5.8S 
•J.  177 

5.1.7 
0. 173 

S.  38 
0.171 

5.14 

O.lbS 

4.80 

O.lfcb 

120.0 

3.81 

0.282 

3.  *33 
0.2b2 

7.4b 

0.247 

b.°2 

0.255 

b.45 

0.224 

b.W 

0.217 

5.b0 

0.208 

5.35 

0.203 

5.0b 

0.137 

4.73 

0.133 

4.55 

0.183 

*4.34 

0.185 

*4.14 
0. 182 

3.'« 

0.178 

3.73 

0.177 

3.b3 

0.174 

1*40.0 

b.09 

0.288 

5.4? 

0.270 

5.30 

0.253 

4.87 

0.244 

4.  bO 
0.234 

4.42 

0.225 

4.13 

0.217 

3-38 

0.212 

3.78 

0.20b 

3.tl 

0.202 

3.45 

0.137 

3.30 

0.183 

3.  lb 
0.181 

3.  CM 
0.133 

2.82 

0.185 

2.81 

0.182 

lbO.O 

u.47 

0.2<W 

4.20 

0.273 

3. 8b 
0.2b4 

3.75 

0.252 

3.5b 

0.242 

3.38 

0.234 

3.22 

0.227 

3.00 

0.220 

2.35 

0.215 

2.32 

0.210 

?.7i 

0.207 

2.bO 

•J.2C3 

2.51 

0.183 

2.41 
0.1  **> 

2.35 

0.183 

2.70 

0.180 

180.0 

3.42 

0.300 

3.24 

0.75*1 

3.C8 

0.271 

2.33 
0.  2bl 

2.80 

0.231 

2.b3 

0.242 

2.5b 
0. 735 

2.4b 

0.223 

7.3b 

0.224 

2.27 

0.213 

2.18 

0.215 

-.11 

0.211 

2.04 

0.207 

1.87 

U.2C4 

l.*» 

0.202 

1.84 

0.188 

200.0 

2.71 

0.30b 

2.53 

0.243 

2.4b 

0.280 

2.3b 
0. 2b3 

2.2b 

0.233 

2.17 
0. 250 

2.08 

o.:-4S 

2.01 

0.233 

1.8*4 
0.  .'33 

1.87 

0.228 

1.80 

0.22*1 

1.74 

0.220 

l.b3 

0.217 

1 .  b*4 
0.214 

1.58 

0.511 

1.54 

0.203 

22*5.0 

2.20 

0.31*4 

2.10 
0. 305 

2.02 

0.287 

1.34 

0.27? 

l.Sb 

0.2-8 

1.80 

0.2bl 

1.75 

0.25? 

l.b? 

0.247 

t.b2 

•5.2*12 

l.Sb 

0.237 

1.51 

0.234 

1.47 

0.230 

1.42 

0.22b 

1.3? 

0.223 

1.34 

0.218 

1.31 

0.2lb 

.'40.0 

1.82 

0.321 

l.T*-* 

0.30? 

l.b3 

0.283 

I.b2 
0. 287 

l.Sb 

0.278 

1.01 
0. 2b3 

1.4b 
0. 2b2 

1.41 

0.23b 

1.37 

0.253 

1.33 

<5.243 

1.23 
0. 2*13 

1.2S 

0.233 

1.72 

0.255 

1.18 

0.231 

1.15 

0.228 

1.12 

0.773 

2*25.0 

1.33 

0.33? 

1.47 

0.313 

1.42 

0.30? 

1.33 

0.23b 

1.33 

0.287 

1.23 

0.273 

1.25 

0.27*1 

1.21 

0.2b3 

1.10 

0.2b2 

1.14 

0.257 

1.11 

0.233 

1.08 

<5.248 

1.(55 
0. 245 

1.03 
0. 241 

1.00 

0.238 

0.87 

0.234 

-SU.0 

1.30 

0.380 

J .  24 
0.323 

1.22 

0.31b 

1.18 

0.30b 

1.15 

0.283 

1.11 
0. 231 

1.03 

0.234 

1.05 
il.  278 

1.02 

0.27? 

0.33 
0.2b  7 

0.3? 
0. 2bS 

0.34 

0.258 

0.32 

0.254 

0.80 
0. 251 

0.88 

0.250 

0.85 

0.24b 

r*J*5.0 

1.12 

0.330 

1.03 

0.333 

1.0b 

0.323 

1.03 

0.313 

1.00 

0.310 

0.87 

0.3*5? 

0.44 

0.2-44 

0.82 

0.288 

0.83 

0.285 

0.37 

0.273 

O.bS 
0. 273 

0.85 
0. 2b3 

0.31 
0.  a>7 

0.73 
0. 2b? 

0.77 
0.  EbO 

0.7b 

0.257 

320.0 

0.86 

O.jh’l 

0.35 

0.332 

0.82 

0.341 

0.85 

0.331 

0.88 

0.322 

0.85 

0.314 

0.83 

O.iCb 

0.81 
0. 300 

0.73 

0.234 

0.77 

0.283 

0.75 

0.237 

0.7*1 

0.283 

1J.72 

0.278 

0.70 
0. 275 

0.b8 

0.271 

0.b7 

0.289 

aExplosive  parameters  for  TNT  were  used  in  these  calculations. 
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Table  B5.  Observed  airblast  overpressures  for  Tugboat  Events  la  and  lb. 


Tugboat  la,  1.0  ton,  Dr  +  =  17.3  ft 

Scaled  depths:  water  =  30  ftyct1^;  coral 

=  143  a  =  0.3425; 

APq  =  6.87  psi 

»./<Dr  +  V 

AP  (psi) 

ap^p0 

23.2 

0.098 

0.0143 

40 

0.042 

0.0061 

102 

0.0155 

0.00226 

132 

0.0096 

0.00140 

198 

0.0078 

0.00114 

234 

0.0040 

0.00058 

371 

0.0032 

0.000466 

423 

0.0027 

0.000393 

479 

0.0018 

0.000262 

Tugboat  lb,  1.0  ton,  Dr  +  Dw  -  17.9  ft 
Scaled  depths:  water  =  31  ft/kt^;  coral 

=  148  tt/kt1/3;  a  =  0.3432, 

APQ  =  6.46  psi 

MDr  *  V 

AP  (psi) 

AP^Pq 

17.3 

0.127 

0.0197 

33.8 

0.041 

0.0063 

93.4 

0.014 

0.00217 

130.0 

0.0137 

0.00212 

187 

0.0082 

0.00127 

230 

0.0064 

0.00099 

364 

0.0038 

0.00059 

467 

0.004  (?) 

Q.00062  (?) 
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Table  B6.  Observed  airblast  overpressures  for  Tugboat  Events  Ic  and  Id. 


Tugboat  Ic,  1.0  ton,  Dr  +  D^.  =  21.7  ft 
Scaled  depths:  water  =  38  ft^t^;  coral 

=  179  ft^kt1^;  a 

=  0.3437;  APQ  =  4.35  psi 

V<Dr  +  V 

AP  (psi) 

*PAP0 

16.4 

0.083 

0.0191 

29.4 

0.040 

0.0092 

78.5 

0.0143 

0.00329 

110.0 

0.0078 

0.00179 

155 

0.0056 

0.00129 

190 

0.0040  (?) 

0.00092  (?) 

298 

0.0030 

0.000690 

340 

0.0028 

0.000644 

384 

0.0036 

0.000823 

Tugboat  Id,  1.0  ton,  Dr  +  Dw  =  25.8  ft 

Scaled  depths:  water  =  45  ftjkt1/9;  coral  =  213  it/kt1^3:  a  =  0.3468;  APQ  =  3.20  psi 


V(Dr  +  D«-> 

AP  (psi) 

AP/APq 

10.0 

0.104 

0.0325 

21.0 

0.038 

0.0119 

62.5 

0.0112 

0.00350 

94.0 

0.0067 

0.00209 

127.0 

0.0042 

0.00131 

163 

0.0029 

0.00091 

255 

0.0026 

0.00081 
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Table  B7.  Observed  airblast  overpressures  for  Tugboat  Event  Ie. 


Tugooat  Ie,  10.0  tons,  Dr  +  D  =  42.9  ft 

Scaled  depths:  water  =  37  ftyct^;  coral  =  163  ffykt*^;  a  =  0.3425;  APQ  =  5.22  psi 


V(Dr  +  V 

AP  (psi) 

AP/APq 

6.5 

0.259 

0.0496 

65.5 

0.0225 

0.00431 

73.9 

0.0185 

0.00354 

103.0 

0.0128 

0.00245 

157 

0.0088 

0.00169 

178 

0.0085 

0.00163 

201 

0.0041 

0.000786 

of  the  peak  water-shock-induced  overpressures  at  intermediate  ranges  (5  to  10  < 

R  /(D  +  D  )  <  100),  and,  if  possible,  photographic  measurements  of  the  peak  water 

W  I*  w 

surface  spall  velocities. 

LIMITATIONS 

MULTIMEDIA  applies  to  single-charge  detonations  in  a  saturated  cratering 
medium  overlain  by  a  water  layer  of  moderate  thickness.  The  shock  transmission 
model,  in  the  cratering  medium,  must  be  known.  At  the  present  time, 

MULTIMEDIA  must  be  considered  a  preliminary  code,  subject  to  revision  and  improve¬ 
ment  of  the  refracted  shock  propagation  model.  Predictions  thus  far  apply  to  a  limited 
range  of  explosive-medium  combinations,  but  this  situation  is  expected  to  improve  as 
more  high-accuracy  shock  transmission  data  for  saturated  media  become  available 
from  SOC-TENSOR  chemical  explosive  modeling  calculations  and  other  sources. 
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I 


FORTRAN 


4 

5 

6 

7 

8 
9 


LIST  OF  MULTIMEDIA  CODE 


1  '*72  AD  G1Q/13  13.9 

OIHENSION  XMED(14)  ,RSUR2(90> , VEL2(90 ) , 0EPW2 (50) ,0EPR2 (50) ,DLP02(50 
1.50) » ALPA2 (50 .50 ) 

INPT =2 
ITPT=3 

DETERMINATION  OF  ALPHA.  WHICH  IS  USED  TP  PERFORM  MONTAN  OR 
THEORETICAL  AIRBLAST  PREDICTIONS.  THIS  PROGRAM  INTEGRATES 


10 

C 

THE  SURFACE  VELOCITY  FIELD  AND  CALCULATES  ALPHA  (AS  WELL  AS 

11 

C 

OTHER  IMPORTANT  QUANTITIES)  FOR  ALL  DETONATIONS  IN  TWO 

-LAYERED 

12 

C 

MEDIA  (SATURATED  ROCK.  ETC.  OVERLAIN  BY  WATER) 

13 

RAP- 1. 0/57.2957  8 

14 

REAO  INPUT  TAPt  INPT ,99, (XMED (N) ,N=1 ,14) 

15 

r9  FORMAT  (13A6.A2) 

16 

WRITE  OU 'PUT  TAPE  ITPT.100 . (XMED (N) , N=l, 14) 

17 

100  FORMAT  (13A6.A2/////) 

18 

REAO  INPUT  TAPE  INPT, 101, RHOCW.CW 

19 

REAO  INPUT  TAPE  INPT ,102, R0CR1,CR1,RTRAN«R0CR2,CR2 

20 

REAO  INPUT  TAPE  INPT ,103, Al, 81 ,RTP, A2, 82 

21 

’.'ll  F  ORMAT  (F5.0.F5.0J 

22 

12  FORMAT  (5F5.0) 

23 

103  FORMAT  (E9. 0, F5. G,F5 .0 ,£9. 0 ,F5.0) 

24 

XNIN=RAD*90.00 

25 

REAO  INPUT  TAPE  INPT ,207, PAMB 

26 

207  FJRHAT  (F8.2) 

27 

00  27  11-1,16 

28 

0£PW=20*II 

29 

0EPW2(II)=0EPW 

30 

IKMAX=II 

31 

DO  27  JJ=3,16 

32 

DEPR=20*JJ 

33 

DEPR2I JJ)=OEPR 

34 

JKMAX=JJ 

35 

WRITE  OUTPUT  TAPE  ITPT.508 

36 

WRITE  OUTPUT  TAPE  ITPT.501 

37 

WRITE  OUTPUT  TAPE  ITPT.502 

38 

500  FORMAT (10H  ANGLE  ,10K  ANGLE  ,10H  SURFACE  ,10H  PRESSURE 

»10H 

39 

1  TRANSMIS  , 10H  PRESSURE  ,10HVERTICAL  ,1 OHCUHULATI VE) 

40 

501  FORMAT (10H  ROCK  ,10H  WATER  ,10H  RANGE  ,10H 

ROCK 

,10H 

41 

1  FACTOP  , 10H  WATER  .10HVELOCITY  .11HVEL  PROFILE) 

42 

502  F ORMAT (10H  OFGREES  ,10H  OEGREES  ,10H  FEET  ,10H 

PSI 

,10H 

43 

1  ,10H  PSI  ,10H  FT/SEC  ,10M  INTEGRAL  ) 

44 

00  216  1=1,86 

45 

X 11=90-1 

46 

ANGLR=RAO*XII 

47 

SR=DEPR/ ( SINF ( ANGLR) 1 

48 

IF  (SR-RTRAN) 217,217,218 

49 

217  RHOCR=ROCRl 

50 

GR=CR1 

51 

GO  TO  219 

52 

218  RH0CR=R0CR2 

53 

CR=CR2 

54 

219  IF(SR-RTR)220,220,22J 

55 

220  PR0C=RH0CR*A1*(SR*»81) 

56 

GO  TO  222 

57 

221  PR0C=RH0CR*A2*(SR»*82) 

58 

222  CSARK=COSF(ANGLR) 

59 

ANGLWsACOSF (CSARK* (CW/CR5) 

60 

SW=0EPW/ (SINF (ANGLW) ) 
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61 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
62 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 
9fc 
99 

100 

ini 

102 

103 

104 

105 

106 

107 

108 

109 

110 
111 
112 

113 

114 

115 

116 

117 

118 

119 

120 
121 
122 

123 

124 

125 

126 

127 

128 

129 

130 


ANGR2=XNIN-ANGLR 
ANGM 2sXNIN-ANGLW 

TANGL=<2.0*RHOCW*COSFfANGR2>)/(RHOCW*COSF<ANGR2)4RHOCR*COSF(ANGW2) 

It 

RSURF=DEPR/ <T  ANF CANGLR) t ♦DEPW/ <TANF ( ANGLW) ) 

A5=0EPR/ ISORTF (RSURFt ) 

A6=SQRTF (TANF  (AMGR2) ) 

A7*l OEPR) ♦ COEPW) * ( IT  ANF  (ANGW2 ) ) *  (1.04CTANF3ANGW2) )  *  IT  ANF  (ANGM2)  >  >  / 
H  <TANC (ANGR2) ) *  1 1.0* (TANF< ANGR2) )•  (T  ANF ( ANGR2 It))) 

A  8*SQRTF  ( A7 1 
A9=A5*A6/A8 

PWAT*TANGL*PROC*  C <A9>  **1.13) 

VEL= (SINF (ANGLW) >*C2 .0  *PWAT/ (RHOCW) > 

RSUR2< I)=RSURF 

V£L2<I)=VEL 

IF  CI-1)  300,300,301 

300  VO=VEL2U> 

ALPHs(V04VEL2<I))*«RSUR2{I)) 

RATIO= VO/ 10 87.0 

GAM= 1.2 

AAA=1. 0 

CCCs-1.0 

B89=GAM*RATI0 

8SQ=GAM*GAM*RATI0*RATI0 

XM*<8B8+S0RTF (BSQ-4. 0* AAA*CCC  >  >/ (2. 0  *AAA ) 

209  AH=XM 

210  OELPO=  CPAMB/ 1013.0) * (0.01841)  *AH*VO 
0LP021XI, JJ1=0ELP0 

GO  TO  379 

301  ALPH=ALPH4<VEl2<mVEL2(I-m  *tABSF(  RSUR2  U1-RSUR2  fl-l)n 
0ELR*RSUR2 (I ) -RSUR2 t  1-1) 

IFCDELRI711,711,379 
711  WRITE  OUTPUT  TAPE  ITPT.714 

714  FORMAT <1CX,82HINV£RSI0N  DUE  TO  BREAK  POINT  APPROX  DELTA  RANGE  IS 
1NEGATIVE  FOR  THE  NEXT  PRINTOUT! 

37  9  X  ANGRsANGLR/R An 
710  XANGWsANGLW/RAO 

WRITE  OUTPUT  TAPE  ITPT ,400,XANGR,XANGW,RSURF,PR0C,  TANGL,PWAT, VEL *A 
1 LPH 

400  FORMAH8F10.3! 

VTEST*VO-50.0*VEL 
IF<VTEST)216, 216.999 

216  CONTINUE 

999  ALPHAS  < ALPM) • <1. 0)/( 2.0*3. 14159* (DEPR40EPW) *VO) 

ALPA2«II,JJ!=ALPHA 

WRITE  OUTPUT  TAPE  IT?T,40i,OEPW,DEPR,AN, VO, DELPO, ALPHA 

401  FORMAT  (/12HWATER  0EPTHs,F8.2,tX,8MFT/KTl/3,3X,12HR0CK  DEPTH*  ,F8. 
12,1X,8HFT/KT1/3,3X,9HMACH  NO  =,F8.5, 3X,17HWATER  VEL  AT  SGZ*,F8.2,1 
1X.7HFT/SEC  /18H0ELTA  P<0>  AT  SGZ=,F9.3,4H  PSI,3X,6HALPHA*,F10.4/// 
1/1 

27  CONTINUE 

WRITE  OUTPUT  TAPE  ITPT, 700 

700  F ORMAT 1 1H1///////5X, 38HSUNMART  OF  DELTA  P(O)  ANO  ALPHA  VALUES1 
WRITE  OUTPUT  TAPE  ITPT, 701 

731  F ORMAT  t///lSX»53HWATER  0EPTHS,FT/KT1 /3,  LISTED  HORIZONTALLY  ACROSS 
1  TOP) 

WRITE  OUTPUT  TAPE  ITPT, 702 

7C2  F ORMAT </2X,39HR0CK  DEPTHS,  FT/KT1/3,LISTE0  VERTICALLY) 

WRITE  OUTPUT  TAPE  ITPT, 703, (OEPW2 (IK ),IK*1, IKNAX) 

703  FORMAT C//16F7.1//I 
00  704  JK=3,JKHAX 

WRITE  OUTPUT  TAPE  ITPT, 705, DEPR2<JK> 

705  F0RMAT</1X, F5.lt 

WRITE  OUTPUT  TAPE  ITPT, 706, <0LP02<IK,JK) .IKsl.IKMAX) 

706  FORMAT (16F7.2) 

WRITE  OUTPUT  TAPE  ITPT, 707, tALPA2 (IK,JK) ,IK*1,IKHAX) 

707  FORMAT C16F7. 3) 

704  CONTINUE 

28  CALL  EXIT 
END 


48 


Appendix  C 

Apparent  Charge  Depth 
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Appendix  C 

Apparent  Charge  Depth 


The  logical  depth  normalization  to  use  in  writing  Eqs.  (13)  and  (16)  would  be  the 
total  apparent  charge  depth  as  seen  from  the  water  layer.  The  total  apparent  charge 
depth,  D  ,  for  a  given  ray  in  the  water  layer  is  the  depth  below  the  water  surface  at 
which  that  ray,  when  projected  back  along  its  direction  of  travel  in  the  water,  inter¬ 
sects  the  vertical  charge  burial  axis.  Due  to  refraction  at  the  rock- water  interface, 
the  apparent  depth,  Dap,  is  not  equal  to  the  true  charge  burial  depth  (Dr  +  Dw).  If 
Cr  or  Ur  >  C^,  the  apparent  depth  will  always  be  greater  than  the  true  depth. 

The  apparent  depth  may  be  determined  using  the  law  of  refraction,  Eq.  (7).  A 
sample  solution  will  be  performed  for  rays  near  the  vertical  direction,  since  the  shock 
pressures  and  spall  velocities  are  highest  in  this  direction.  From  the  law  of  refraction, 

cr 

cos  er  =  cos  ew  * 


(In  the  inelastic  region,  Ur  should  be  substituted  for  Cr  in  this  and  subsequent  equa¬ 
tions.) 

For  rays  near  the  vertical  direction. 


tan  ew  =  tan  Q 


rTT 
w 


The  total  apparent  depth  D  for  a  ray  (near  the  vertical  direction)  that  reaches  the 

ap 

water  surface  at  surface  range  is  then  given  by: 


Dap  ■  ta"  W 


"  torn  8r<Cr/Cw?< R„> 


Substituting  from  Eq.  (10)  for  Rw  gives: 


=  *  D»  (®n£ 


=  constant. 


Let  us  consider  only  the  portion  of  the  far-field  overpressure  pulse  produced 
by  these  rays  near  the  vertical  direction.  Using  the  above  result  (Dap  =  constant)  in 
Eqs.  13  and  IS, 


a 


V( water,  at  Rw)  dRw  , 
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where  ARw  is  a  range  close  to  SGZ  (i.e.,  AR^  —  0);  and 


ap“apo«-r E  =apo 


0  ap 


V( water,  at  Rw)  dR 


.)(£*)  • 
/  '  w  / 


The  D  terms  in  this  equation  for  AP  cancel,  allowing  the  use  of  any  desired  depth 
ap 

normalization  in  place  of  Dap,  including  the  true  depth  Dr  +  Dw  (*hich  happens  to  be 
the  most  convenient).  Note  that  the  same  depth  normalization  must  be  used  in  both  the 
equations  for  a  and  AP.  More  generally,  it  can  be  shown  that  this  result  holds  for  that 
part  of  a  far-field  overpressure  pulse  AP  produced  by  any  small  segment  of  the  surface 
velocity  field,  whether  or  not  the  ray  producing  that  segment  is  near  the  vertical 
direction.  Thus,  the  result  applies  for  the  entire  integrated  surface  velocity  field  and 
for  a,  above. 
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